^ r 


/K /-</3 

% I 3 A. 


Final Technical Report 


to 


National Aeronautics and Space Administration 


NASA Grant No. NAG 5-919 


Remote Measurement of Soil Moisture Over Vegetation 
Using Infrared Temperature Measurements 


Submitted by: 

Dr. Toby N. Carlson 

Department of Meteorology 
The Pennsylvania State University 
University Park, PA 16802 


January 1991 

ASA-CR-ld73*2) REMOTE MEASUREMENT Of SOIL 
I5TURF OVER VEGETATION USING INFRARED 
TEMPERATURE MEASUREMENT* Final Technical 
Report (Pennsylvania State Univ.) i 32 p 

CSCL 02F Gi/43 


N91-I561 4 

Unci as 
0326318 


1. Introduction 


This final report summarizes research work performed between the inception of 
this grant in 1986 and the present. The underlying motivation for all the work was to 
develop better methods for remote sensing of surface evapotranspiration, soil moisture 
and fractional vegetation cover. Details and results are to be found in several journal 
articles and four theses. In the interest of bre vity, we will refer directly to these 
publications, which are attached as appendices, and to their figures. These papers are, 
respectively: Appendix A (Carlson a nd Buff um, 1989); Appendix B (Lynn and Carlson, 
1990); Appendix C (Carlson et al., 1990); Appendix D (Carlson et al., 1991). Additional 
publications generated by this grant and referred to in subsequent sections consist of 
reviewed publications (Carlson and Lynn, 1991; Carlson, 1991; Andre et al., 1988), one 
publication under review (Carlson and Lynn, 1991) and theses (Buffum, 1988; Lynn, 
1989; Belles, 1990). 

Objectives of our research were (1) to further develop a model of water 
movement through the soil\plant\atmosphere system, (2) to use this model, in 
conjunction with measurements of infrared surface temperature and vegetation index, to 
calculate soil water content and vegetation fraction, (3) to determine the magnitude of 
radiometric temperature response to water stress in vegetation, (4) to show at what 
point one can detect that sensitivity to water stress, and (5) to determine the practical 
limi ts of the methods. A new approach, undertaken during this past year, was to 
develop a hydrological model that can be used to calculate soil water content versus 
depth given conventional meteorological records and observations of vegetation cover. 
We will now outline the results of these initiatives. 


2. Projects Undertaken 

(a) Large-scale evapotranspiration measurements 

We investigated the method developed by Jackson et al. (1977) to determine daily 
evapotranspiration from a single radiometric surface temperature made near mid day. 
The idea was to use our plant/soil/atmosphere model to simulate radiometric surface 
temperatures and to see if we could reproduce the experimental results and find the 
same numerical values for the coefficients in the formulae given in the literature 
(Jackson, 1977; Seguin and Itier, 1983). Our results (summarized in Appendix A 
(Carlson and Buffum, 1989) and in Buffum (1988) show that the model yields very 
similar coefficients to that of Jackson et al. (1977) for vegetation and Seguin and Itier 
(1983) for bare soil. Surface roughness and wind speed also exert an important effect 
on the coefficients. 

We show, however, that it is not necessary to specify more than two or three wind 
speed and roughness categories in order to achieve a useful accuracy. The method 
confirms that it is theoretically possible to determine the net soil water loss by 
evapotranspiration using a midday radiometric surface temperature. However, since the 
coefficients for vegetation and bare soil are so different, one should also have knowledge 
of the vegetation fraction. Carlson and Buffum (1989) also present a variation of the 
method which does not require air temperature. In it, one substitutes for the air 



temperature a second surface radiometric temperature during the morning to obtain a 
rate of morning increase in temperature. 

b) Modeling the plant canopy 

Much of our modeling effort has been devoted to improving the vegetation 
component of the model. A serious deficiency in using radiometric temperature to infer 
soil water content in the presence of vegetation is that, according to both observations 
and simulations, the radiometric surface temperature of a vegetation canopy tends to 
vary only slightly as a function of soil water content over a wide range of the latter. 
Typically, vegetation canopies appear "cool" (within a couple degrees centigrade of air 
temperature) until some low value of soil water content is reached, below which the 
leaves may warm rapidly with decreasing soil water content. This suggests, that plants 
may transpire at nearly an optimum rate until urgent water stress signals are received 
from the roots. 

In order to approach this problem it was necessary to rethink our parameterization 
of stomatal resistance. Many stomatal resistance models, such as the one we had 
previously used, either fail to account for the lack of sensitivity of stomatal behavior to 
soil water content or are unable to describe the effects of transient stress in which leaf 
temperature can temporarily rise in response to excessive atmospheric demand in the 
face of reduced soil moisture. 

Our present plant model incorporates a more mechanistic approach than before. 

By mechanistic we mean that the model accounts for internal plant mechanisms, rather 
than simply calculating the surface energy fluxes, as in some of the bulk 
evapotranspiration methods. At the same time we have made the model as 
deterministic as possible. By deterministic, we mean a mathematical description of the 
evapotranspiration as a function of external conditions or internally calculable or 
prescribed variables. To some extent, the idea of a mechanistic model is an elusive one, 
since any model must inevitably simplify the infinitely complex biochemical processes 
that constitute the plant mechanisms. Nevertheless, there is at present a hierarchy of 
mechanistic complexity in plant models, extending from simple turbulent transfer 
equations that take no account of the water supply in plants to ones which treat the 
biological functions of the leaves, stem and roots. 

The new stomatal resistance model reported by Lynn and Carlson (1990; Appendix 
B) simulates a rich variety of plant behavior, including transient effects of water stress 
caused by a excess atmospheric demand. It does this by allowing stomatal resistance to 
depend on the (epidermal) leaf water potential, subject to a threshold leaf water 
potential below which the stomatal resistance increases rapidly with decreasing leaf 
water potential. Although the dependence of stomatal resistance on leaf water potential 
is questioned by many plant scientists, there is considerable experimental evidence to 
suggest that the loss of turgor pressure in the leaf epidermis (guard cells) causes a rapid 
increase in stomatal resistance with decreasing epidermal water potential below a 
threshold epidermal water potential and a corresponding decrease in transpiration on 
the scale of a field (Federer, 1980; Nizinski and Saugier, 1989; Kateiji, 1990; Carlson et 
aL, 1991; Appendix C). Our simulations (Figure 1) show that the epidermal water 
potential and vapor pressure deficit effects on stomatal resistance are intertwined. 



Inclusion of the threshold leaf water potential in the stomatal resistance model 
allows the latter to simulate transient periods of plant water stress, which can occur even 
when the soil is well watered provided that the atmospheric demand is large. Figure 14 
in Appendix B shows two evapotranspiration curves made over adjacent corn fields, one 
irrigated and the other non-irrigated. Neither soil was excessively dry. Both 
evapotranspiration curves are identical except during a few hours when the 
evapotranspiration over the dry field exhibited a plateau (or"hat M ), signifying a limitation 
in the plant’s rate of water loss. The hat can be simulated by imposing a threshold 
epidermal water potential (^J, in this case of -13 bars. Figure 9 in Appendix B shows 
the sensitivity of stomatal resistance and transpiration to the choice of threshold leaf 
water potential. (Henceforth, unless otherwise stated, we will interchange the terms, 
evapotranspiration and transpiration, realizing that the two are similar over dense 
vegetation.) Not surprisingly, the threshold leaf water potential is a crucial parameter in 
the model. Typically, this value of ^ is about -15 bars in most temperate crops, but it 
can vary from -9 bars to -25 bars, depending on species. 

The transpiration plateau occurs as follows: a decrease of the epidermal leaf 
water potential (to more negative values) occurs during the morning because the 
transpiration is increasing. Depending on atmospheric demand (essentially solar flux 
and vapor pressure deficit) and the soil water content, the epidermal leaf water potential 
may reach the threshold value. The greater the atmospheric demand, the larger the 
threshold water potential or the lower the soil moisture the more likely it is that the leaf 
water potential will reach the threshold and the longer the threshold will last 

When the epidermal water potential reaches the threshold, stomatal resistance 
increases rapidly with a further decrease in leaf water potential (reversing a decrease 
with time); this may lead to a "bump" in the stomatal resistance, as shown by both the 
(HAPEX) measurements and the simulations in Figure 12 of Appendix B. Because of a 
concomitant increase in vapor pressure deficit (between the leaf and its surroundings), 
due to the warming of the leaves, the increase in stomatal resistance is offset by an 
increase in atmospheric demand; the resultant effect of reaching the threshold leaf water 

potential on transpiration is a plateau rather than a collapse of the transpiration or even 
a decrease; the complex feedback between various mechanisms is captured by the 
model, showing that the plant makes optimum use of water in order to avoid becoming 
desiccated with increasing atmospheric demand. 

The point to be made here is that water stress is relative, and the rise in leaf (and 
therefore canopy) temperature depends as much on atmospheric demand in the face of 
plant constraints as on the soil water content. Before complete wilting occurs, however, 
the transient stress, manifested by the "hat" results in an augmentation in the increase of 
leaf temperature, thereby affording an early warning of plant water stress through its 
detection by radiometer; this subject is addressed in the paper by Carlson et al. (1991; 
Appendix D). 

A unique feature of the stomatal resistance model is the way in which vapor 
pressure deficit is specified. Without going into detail, we refer the reader to Eqn. (9) 
in Appendix B. Briefly, this equation relates the vapor pressure deficit to the difference 


in water potential between the mesophyll and the epidermis. This formalism allows the 
epidermal and mesophyll leaf water potentials to vary in the opposite sense, as shown in 
Figure 1. It also allows transpiration to begin decreasing with increasing vapor pressure 
deficit as the latter becomes large; this is known as the feedforward effect (Farquhar, 
1978). Inclusion of a large vapor pressure deficit sensitivity, as in the case of soybeans, 
results in the evapotranspiration plateau becoming concave. This is shown in the results 
of simulations and measurements of evapotranspiration made at Penn State (Figure 2). 
Further details are contained in a thesis by Lynn (1989). 

i) Stomatal resistance workshop 

A workshop conceived and organized by the PI on the subject of stomatal 
resistance was held at Penn State, 10 - 13 April, 1989. Partially funded by this grant, the 
workshop hosted about 35 invitees for the purpose of discussing methodology and need 
for improved measurement and modeling of stomatal resistance. Besides the intensive 
and lively discussions, the workshop resulted in a set of proceedings (available on 
request) and a collection of 13 reviewed articles to be published during 1991 in 
Agricultural and Forest Meteorology. While not contributing directly to these articles, 
the PI served as acting editor and has written a preface to the volume. 

c) Partial vegetation cover 

Spatial variations in surface radiometric temperature may be due to spatial 
variations in vegetation cover, rather than in soil water content The reason is that 
upwelling thermal radiance is emitted from an arbitrary mix of warm, dry soil and cooler 
vegetation. Thus, remote radiometric surface temperature measurements can yield 
misleading information on the horizontal distribution of soil water content, depending on 
the amount of bare soil visible to the radiometer, because the vertical variation of soil 
water content becomes artificially mapped onto the horizontal. 

The problem is to devolve both the vegetation fraction and the horizontal soil 
water content from remote measurements. To this end we have made use of radiance 
measurements at solar wavelengths, specifically the normalized difference vegetation 
index (NDVI). We were obliged to spend a considerable amount of time developing a 
partial vegetation canopy component in the boundary layer model, where bare soil and 
vegetation modes meld within a framework of a common mixing layer and substrate. 



Our method of analysis is as follows: we first determine asymptotic surface 
temperatures (the extrema) for sunlit bare soil and sunlit vegetation by plotting the 
NDVI versus the radiometric surface temperature. In many cases this variation shows a 
narrow band of data points extending from high NDVI and low radiometric surface 
temperature to low NDVI (close to zero) and high surface temperature, as shown in 
Figure 6 of Appendix C. We call this narrow strip of data points the "axis of variation". 
The end points of the axis of variation correspond approximately to values for sunlit 
bare soil and a full, sunlit vegetation cover. Note that the variation of NDVI versus 
radiometric surface temperature flattens at high NDVI, signifying that at some point the 
variation ceases to increase significantly with increasing vegetation amount. Carlson et 
al. (1990; Appendix C) discuss the fact that NDVI increases with increasing vegetation 
amount, but effectively saturates (increases insignificantly) with increasing vegetation 
amount above a leaf area index of about 3 for many vegetation types. 

We next assume that the lower temperature extremum, where NDVI no longer 
significantly increases with decreasing surface temperature, corresponds to a leaf area 
index of 3 and to 100% vegetation cover. The high temperature extremum corresponds 
to sunlit bare soil. Accordingly, we deduce the bare soil (surface) water content and the 
vegetation (root zone) water content with the aid of the aforementioned boundary layer 
model, by forcing the simulated temperatures to match measured temperature extrema 
in two simulations; these extrema correspond to that for bare soil and that for full 
vegetation cover (leaf area index = 3), for which we derive the root zone water content. 

The next step in the analysis is to simulate the radiometric surface temperature in 
which we vary vegetation fraction but keep the substrate water contents (for surface and 
root zone) equal to the values derived in the previous step. The result is an estimate of 
surface radiometric temperature versus vegetation fraction, from which we derive the 
relationship between NDVI and vegetation fraction from the axis of variation. In 
principle, we can then produce a map of vegetation fraction and the surface energy 
fluxes. 

Although there is some artistry and guesswork involved, the method offers some 
prospect for remotely measuring soil water content at two levels (surface and root zone), 
vegetation fraction and surface energy fluxes over sparse or partial vegetation cover. 

The method is currently being tested by our co-worker, Rob Gillies, using AVHRR 
satellite images of an urban area in Britain (see below). 

d) Radiometric detection of transient water stress in vegetation 

A question posed in our most recent paper (Carlson et al., 1991; Appendix D) is 
whether transient water stress associated with the transpiration plateau can be detected 
radiometrically. Since the onset of the hat signifies a period of plant water stress and 
may presage a later collapse of transpiration with further soil drying, detection of its 
effect on radiative surface temperature is of some interest We began this project by 
using micrometeorological measurements (surface energy fluxes, radiometric 
temperature, soil moisture) taken during two growing seasons (one in corn and the other 
in soybeans) at Rock Springs, a Penn State University agricultural site. We confined our 
attention to days in which there was (1) dense vegetation (evapotranspiration 
approximately equals transpiration), (2) clear skies, and (3) functioning instruments. In 



the end, we found only one period meeting these criteria that lasted more than one day; 
the measurements pertain to com during July 1987. 

Evapotranspiration measurements during this period show a plateau which 
lengthened from about two hours on the first day to about 5 hours after 3 days of soil 
drying. A comparison of measurements and simulations are shown in Figure 1 of 
Appendix D. The model was capable of simulating the transpiration plateau, including 
its increase in duration and the decrease in maximum evapotranspiration during the 
period; the threshold water potential in the model was set at -15 bars. 

Having successfully simulated the behavior of the hat, we proceeded to make 
additional simulations that suggest how the hat may have evolved with time. We wanted 
to show how well the rise in canopy temperature could be detected radiometrically. 
Figure 9a (Appendix D) shows the simulated change in radiometric surface minus air 
temperature differences at 1300 local time as a function of plateau length and date in 
July 1987 (numbers above curve). Let us assume that the basic uncertainty in measuring 
a change in surface minus air temperature due to soil water content at a single point is 
about two decrees C. Since the soil was almost at field capacity on July 4th (the last 
day with significant rain before the end of the month), the departure of the canopy 
radiometric temperature from that for a soil at field capacity would not have exceeded 
two degrees C until about July 22. Therefore, the period of stress marked by the 
plateau would not have been detectable until a period of more than two weeks of drying 
had passed. By that time the soil water content would have decreased to less than 50% 
of field capacity and the length of the plateau exceeded 3 hours. Interestingly, however, 
simulations suggest that the plateau begins to form after just a few days when the soil 
water content was still rather high, about 70% of field capacity. This early evidence of 
water stress occurs because of high atmospheric demand rather than low soil moisture. 

Another interesting conclusion suggested by these simulations is that the plateau 
would not occur if the atmospheric forcing were decreased by the reduction of sunlight 
associated with cloud cover (Figure 13 in Appendix D). When the net radiation was 
diminished to that of a broken cloud cover, the plateau vanished completely. Moreover, 
we also found that the plateau was effectively masked when the fractional vegetation 
cover was reduced to about 50% (Figure 11a in Appendix D). This suggests that 
detection of plant water stress using radiometric temperatures is dependent on a premise 
of clear skies and dense vegetation; variations with time in cloudiness or in space of 
vegetation fraction can completely obscure the water stress signal. In such cases of low 
atmospheric demand, the radiometric surface temperature is more likely responding to 
variations in surface (rather than root zone) soil water content. 

The results of this study are very intriguing. They suggest, as has Katerji (1990), 
that the threshold leaf water potential can be inferred from evapotranspiration 
measurements by observing the plateau and then simulating the length of the plateau 
and the evapotranspiration level during the plateau. Further details of the field 
experiment are contained in a thesis by Belles (1990). 



e) The effects of plant capacitance 


A paper by Carlson and Lynn (1991) discusses the effect of including plant 
capacitance (water storage) in the model reported by Lynn and Carlson (1990). We 
show that the transpiration plateau can be mitigated or even removed by the effects of 
capacitance. Surprisingly, our results show that the effects of capacitance on either 
transpiration or stomatal resistance is negligible except when there would otherwise be a 
transpiration plateau. We also show that capacitance may constitute an important effect 
in the field- scale transpiration, but probably only for large vegetation (e.g. trees). Thus, 
trees and other large plants may have the capability of eliminating the transpiration 
plateau as the result of internal water storage. This paper is currently under review. A 
brief discussion of the this paper and those contained in appendices A, B and C are 
contained in a short review paper by Carlson (1991). 

f) Other projects 

i) A "bottom up" method for calculating soil water content 

As we have discussed above, radiometric temperature measurements over vegetation 
can not detect deep soil drying until the soil water content has decreased below about 
50% of field capacity. Nevertheless, variations in soil water content between 50 and 
100% of field capacity are important for water budget analyses and for initializing 
regional and climate prediction models. Because of limitations in the radiometric 
method, we have taken a complimentary approach, which we refer to as the "bottom up" 
method, as opposed to the radiometric approach, which we refer to as the "top down" 
method. 

The bottom up method uses routine meteorological and surface observations in 
conjunction with a simple hydrological model to calculate soil water content. The model 
takes precipitation, cloud cover, vegetation height and fractional cover, maximum and 
minimum temperature, the depth of the water table and some knowledge of the soil type 
and it calculates evapotranspiration, runoff, diffusion of water in the soil, and infiltration 
in the process of determining the vertical profile of water content in the soil. This 
model will be reported in a forthcoming MS thesis by William Capehart. 

Figure 3 shows some preliminary results with this model, specifically the surface 
volumetric water content simulated for the same Penn State agricultural site referred to 
in Carlson et al. (1991). We wanted to find out how long it took for the simulations to 
"forget" the initial conditions, which are usually not known. We started the simulations 
with three differing soil water content profiles on 1 January 1987, By July, 1987, the 0 - 
10 cm values have merged, indicating that the model is no longer sensitive to the choice 
of initial conditions after several months. The significance of these results is that one can 
simulate the substrate water contents by starting the model with any reasonable 
substrate water profile several months before the target date. The bottom up method is 
being tested against measurements of soil water content using data collected via gypsum 
blocks during the 1987-88 Penn State field program and with measurements made dining 
the recent MAC Hydro experiment (referred to below). (At present we are conducting 
sensitivity tests to investigate the impact of initial conditions on the results, as well as 



the effects of the water table level, the soil type, and the vegetation height and 
fractional vegetation cover on the vertical profile of soil water content.) 

ii) Testing the partial vegetation cover analysis 

We are currently analyzing images of radiometric surface temperature and NDVI 
over Newcastle, England, in order to derive patterns of vegetation fraction and surface 
energy fluxes. This project will test our method of estimating partial vegetation fraction, 
in this case over a region where vegetation fraction and therefore the evapotranspiration 
has been modified substantially by urbanization. The technique to be applied is 
identical to that referred to above in section (c). Results using AVHRR imagery are 
remarkably similar to those for HAPEX, which pertain to a small-scale agricultural area. 
Figure 4 shows the variation of radiometric surface temperature versus NDVI for the 
Newcastle area on 21 May, 1989. Note that the configuration of data points closely 
resembles that in Figure 6 in Appendix C. The warm temperature extremum, about 
34X1 pertains to a pure urban surface and to a very low surface soil water content. At 
the cool end, the variation bends to the left below about 23°C, approximately where the 
vegetation fraction reaches 100%. Accordingly, we can map the vegetation fraction as a 
function of NDVI between 0 and 100% percent, assuming that the vegetated fraction up 
to 100% corresponds to a leaf area index of 3.0. 

iii) Mahantango field program ( MAC Hydro experiment) 

The Mahantango Watershed field program was conducted by the Agricultural 
Research Service, USDA (Beltsville, MD) and by NASA at Goddard Space Flight 
Center for the purpose of evaluating the usefulness of microwave remote measurements 
of soil water content and the modeling of hydrologic and biogeochemical processes. 
Mahantango Creek watershed, located in central Pennsylvania, is an instrumented 
research basin managed by the USDA, Agricultural Research Service, Northeast 
Watershed Research Center at University Park, PA. We are using the remote 
measurements of surface radiometric temperature and NDVI taken from aircraft and the 
in-situ measurements of vegetation height and fraction and surface meteorological data 
in order to calculate soil water content from the bottom-up and top-down methods 
referred to above. Specifically, we use the remote measurements and our hydrology 
model to calculate soil water content and vegetation fraction. These results will be 
compared with in situ measurements of soil water content and vegetation fraction and 
with soil water content measurements obtained from microwave data; the latter are 
currently being analyzed by the USDA. Preliminary data processing has only just begun; 
further efforts will be made in conjunction with scientists at USDA and at Penn State. 



3. Synthesis and future work 

Sensitivity tests with our vegetation model suggest that water stress in vegetation 
may be detectable for a period of a few days prior to wilting, provided that the 
atmospheric demand is large, an appropriate threshold leaf water potential exists and 
plant storage capacity is not large. What is commonly measured by radiometer are 
variations in the soil surface temperature, depending on the density of the vegetation 
canopy. On the other hand, horizontal variations in radiometric surface temperature 
over sparse vegetation may simply reflect horizontal variations in canopy density rather 
than horizontal variations in soil water content This variation can be used to map the 
vegetation fraction, rather than the soil water content. 

In the figure, we hope to clarify the issue of radiometric measurement of soil 
drying and to show how it can be applied over sparse vegetation. To do this it is 
necessarq to spend some time testing our present methods using available field 
measurements. The MAC-Hydro experiment offers the best opportunity to do so, using 
data obtained over a typical watershed. We will also analyze vegetation fraction for 
Newcastle, England. The latter analyses will be compared with available land use 
information. Modeling efforts will be confined to making the boundary layer-plant 
model user friendly and to include carbon dioxide fluxes in the plant component. 
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Figure 1. Simulated mesophyll leaf water potential (V>^) (bar; crosses) and 

epidermal leaf water potential (V> e ) (plus signs) versus time for soybean 

simulations presented in Figure 2. Vertical arrows denote duration of 
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threshold epidermal leaf water potential. 
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Figure 2. Simulated transpiration (L E ■ solid line labeled in Wm ) and 
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evapotranspiration (L^E; crosses) and measured evapotranspiration (MEASURED; 
triangles) versus time for soybeans on 6 July, 1988, at the Penn State 
agricultural field site. Vertical arrows denote the duration of the 
evapotranspiration plateau in the simulations. 

Figure 3. Simulated surface (0 - 10 cm) volumetric water content versus date 
(1987) using meteorological observations taken at Rock Springs for three 
differing initial soil water profiles, 20, 44 and 60% of field capacity (0.34 
by volume) . 

Figure 4. Variation of NDVI versus radiometric surface temperature (°C) for 
area around Newcastle, England, 21 May, 1989, from AVHRR of NOAA-11 satellite; 
the footprint for each data point is about 1.1 km. 
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Figure 1. Simulated mesophyll leaf water potential (^.) (bar; crosses) and 
epidermal leaf water potential (ij> ) (plus signs) versus time for soybean 
simulations presented in Figure 2. Vertical arrows denote duration of 
evapotranspiration plateau. Horizontal dashed line labeled ^ denotes the 
threshold epidermal leaf water potential. c 
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Figure 2. Simulated transpiration (L E _; solid line labeled in Wm" 2 ) and 
evapotranspiration (L E; crosses) and measured evapotranspiration (MEASURED; 
triangles) versus time for soybeans on 6 July, 1988, at the Penn State 
agricultural field site. Vertical arrows denote the duration of the 
evapotranspiration plateau in the simulations. 
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Figure 3. Simulated surface (0 - 10 cm) volumetric water content versus date 
(1987) using meteorological observations taken at Rock Springs for three 
differing initial soil water profiles, 20, 44 and 60% of field capacity (0.34 
by volume) . 
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Figure 4. Variation of NDVI versus radiometric surface temperature (°C) for 
area around Newcastle, England, 21 May, 1989, from AVHRR of NOAA-11 satellite; 
the footprint for each data point is about 1.1 km. 
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APPENDIX A 



REMOTE SENS. ENVIRON. 29:197-207 (1989) 


On Estimating Total Daily 
Evapotranspiration from Remote Surface 
Temperature Measurements 


Toby N. Carlson 

Department of Meteorology, Pennsylvania State University 

Martha J. Buffum 

Computer Sciences Corp., Falls Church, Virginia 


Al method for calculating daily evapotranspira- 
tion from the daily surface energy budget using 
remotely sensed surface temperature and several 
meteorological variables is presented. Values of the 
coefficients are determined from simulations with 
a one-dimensional boundary layer model with veg- 
etation cover. Model constants are obtained for 
vegetation and bare soil at two air temperature 
and wind speed levels over a range of surface 
roughness and wind speeds. A different means of 
estimating the daily evapotranspiration based on 
the time rate of increase of surface temperature 
during the morning is also considered. Both the 
equations using our model-derived constants and 
field measurements are evaluated, and a discussion 
of sources of error in the use of the formulation is 
given. 


INTRODUCTION 

A desirable goal in remote sensing is to determine 
the total amount of water evaporated from a sur- 
face during 1 day. Use of remotely sensed temper- 

Address correspondence to Dr. Toby N. Carlson, Dept, 
of Meteorol., Pennsylvania State Univ., 503 Walker Bldg., 
University Park, PA 16802. 
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ature for inferring surface evapotranspiration is 
discussed by Bartholic et al. (1972), Jackson et al 
(1977), Soer (1980), and Taconet et al. (1986) for 
regions the size of a single field and by Price 
(1982), Carlson (1986), Flores and Carlson (1987), 
and Wetzel and Woodward (1987) for larger-scale 
areas. These models are generally based on the 
surface energy budget in combination with gov- 
erning equations for the fluxes of water, sensible 
heat, and radiative energy between the surface 
and atmosphere and in the substrate. 

The problem with complex models is that they 
require a detailed set of initial conditions for the 
atmosphere, surface, and substrate in order to 
obtain a solution. As a result, semiempirical meth- 
ods have been proposed for estimating the total 
daily evapotranspiration integrated over a 24-h 
period (LE^) from remotely sensed surface tem- 
perature measurements. One such formula, first 
proposed by Jackson et al. (1977), is 

L£ 24 = R n24 = A + B(r 0 -r o ), (1) 

where R n24 is the net radiation integrated over a 
24-h period (in units of cm), A is a constant, T 0 is 
the radiometric surface temperature, and T a is the 
air temperature. Both T 0 and T a are measured 
near solar noon. The slope B and intercept A are 
determined by linear least squares fit to their data. 
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Choice of measurement levels is not always clearly 
noted in the literature. Seguin and Itier (1983) 
chose 2 m as a standard height and Jackson et al. 
(1977) a level of 1.5 m for T a . Conventional meteo- 
rological observations customarily assign the level 
of surface wind speed measurement at 10 m; in 
micrometeorologica] operations, wind speed is of- 
ten measured closer to the ground. 

An advantage of (1) is its simplicity, requiring 
minimal amounts of ground-based meteorological 
data. More specifically, only three variables (fi n24 , 
T (l , and T 0 ) need to be measured. Thus, the for- 
mula is convenient to use in conjunction with 
satellite measurements of radiometric surface tem- 
perature, such as made by the geostationary earth 
satellite (GOES) or by the NOAA AVHRR at 
2 a.m. and 2 p.m. overpass times (Lagouarde, 
1988). Since 24-h integration of sensible heat stored 
in the substrate is likely to be close to zero, (1) 
expresses the 24-h integrated surface sensible heat 
flux into the atmosphere. 

The B-method, as we will henceforth refer to 
the use of equations such as (1), arose from analy- 
sis of a daily water stress index, which was defined 
by Jackson et al. (1977) for surface and air temper- 
atures made near 1300 local standard time (LST). 
Jackson et al. found that a plot of surface-air 
temperature differences over irrigated wheat ver- 
sus H n24 — LE 24 yielded a straight line with B = 
0.064 cm°C" 1 and A = 0. Seguin et al, (1982) 
collected data over large homogeneous areas in 
southeastern France and found (1) required B = 
0.025 and A = 0.1. 

Seguin and Itier (1983), using the results of 
Itier and Riou (1982), obtained a modified version 
of (1), which is 

R n2A -LE^B{T 0 -T a )\ (2) 

where n is an exponent. Seguin and Itier found 
that n was equal to 1.0 for stable or near neutral 
conditions and 1.5 when the stratification was 
unstable. They also showed that the value of B 
depends upon the roughness length z 0 and the 
wind speed. Consequently, Seguin and Itier (1983) 
suggested that the imposition of a single value for 
B and n, even for various climates and types of 
surfaces, may be unacceptable, although specific 
values of B may be chosen according to certain 
broad categories of roughness, static stability, and 
vegetation type. They concluded that, for “medium 
rough” surfaces, it was worth assigning only two 


values of B, one for stable and the other for 
unstable conditions. They suggested a value of 
B = 0.025 and n - 1 in (2) for unstable conditions. 

An alternate form of (2) was suggested by 
Nieuwenhuis et al. (1985). They eliminated R n24 
and T a and substituted the potential evapotranspi- 
ration and a corresponding reference canopy tem- 
perature. Their modified form of (2) depended on 
a single measurement of temperature with respect 
to a modeled temperature for a potentially trans- 
piring crop. L£ 24 was calculated as a difference 
from a modeled daily integrated potential evapo- 
transpiration, which was determined for a com- 
pletely wet surface using a boundary layer model 
[the TERGRA model of Soer (1980)]. Their equa- 
tion is otherwise similar to (2), except that it is 
expressed in terms of a ratio of actual to potential 
evapotranspiration. This ratio is similar to the 
moisture availability parameter defined by Garlson 
(1986), except that the latter pertains to an instan- 
taneous ratio of evapotranspiration to potential 
evapotranspiration . 

Nieuwenhuis et al. (1985) demonstrate the ef- 
ficacy of the B-method for determining LE ^ over 
individual fields; they also showed that the values 
of the coefficients are sensitive to surface rough- 
ness and meteorological conditions. More impor- 
tantly, Nieuwenhuis et al. (1985) were the first to 
use a full boundary layer model to generate the 
coefficient B. Rambal et al. (1985) performed a 
study similar to ours in which values of B were 
determined from a boundary layer model with 
vegetation component for the case of n = 1.5. They 
found that B increases with increasing height and 
density of tree cover. More recently, Lagouarde 
and Brunet (1988) investigated the variation of A 
and B in (1) to changes in surface roughness using 
a planetary boundary layer model. They show that 
B varies rapidly with changing surface roughness 
when the latter is in the range between 0.2 and 10 
cm, particularly when the surface -air temperature 
differences are large. 

We have also used a sophisticated planetary 
boundary layer model to investigate the response 
of the coefficient B and the exponent n in (2) to 
changes in surface roughness, wind speed, and 
vegetation. This model (henceforth referred to as 
CM) has been described by Carlson (1986) and 
Taconet et al. (1986). Until now, the primary 
practical application of the CM has been in remote 
sensing of the surface energy fluxes; in this con- 
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text, the model is inverted in conjunction with 
remote measurements of surface temperature and 
other parameters to obtain a solution for the sur- 
face energy fluxes and soil moisture availability. 

A different form of (2) which does not require 
measurement of an air temperature for calculating 
the sensible heat flux term is proposed here. 
Wetzel et al. (1984) show that the soil moisture 
(and therefore the evapotranspiration) is most sen- 
sitive to the rate of temperature rise during the 
morning (e.g M between 8 and 10 local time). As a 
result the following equation will be examined: 

R n2i -LE^B'(AT/M) n ' (cm day' 1 ). (3) 

Here (AT/ At) is an average rate of temperature 
rise during the morning (expressed in °C per hour) 
and B' and n' are constants. It should be empha- 
sized here that (3) represents an intuitive leap 
from (2). Justification of (3) is based on data 
generated by the CM, which shows that the rate of 
surface temperature change with time is highly 
correlated with the integrated daily sensible heat 
flux, as expressed by (I) or (2). In omitting refer- 
ence to T a , (3) is suggestive of the equation tested 
by Nieuwenhuis et al. (1985). 

The purpose of this paper is to examine the 
B-method with regard to simulations of the param- 
eters in (2) and (3) by the CM. We will show that 
the parameters B and n in these equations depend 
on surface roughness, wind speed, reference height, 
and vegetation. Two sets of simulations were made 
of B and n (and B' and n' ): at 2 and 6.4 m, 
respectively, for temperature and wind speed and 
at 50 m for both temperature and wind speed. 
Choice of the 2 and 6.4 m levels is dictated by the 
availability of field measurements for testing the 
B-method. Our object is not to add to the bewil- 
dering profusion of values for B and n currently 
available in the literature but to illustrate the sensi- 
tivity of B and n to ambient and surface influ- 
ences and to suggest an alternate form of the 
B-method. Therefore, a part of this paper is de- 
voted to demonstrating the validity of (3). 

METHODS 

The CM describes the water and energy exchanges 
between the substrate, the layer of vegetation, an 
atmospheric surface layer and a mixing layer. Sur- 
face layer similarly laws are employed to compute 
the surface energy fluxes with corrections for 


changes in static stability. An account is taken of 
differences between the thermal roughness length 
for heat and for momentum. Results are deter- 
mined by a complex nonlinear interaction with 
evolving static stability and ambient conditions in 
a multilayered atmosphere and substrate. 

Numerous simulations performed with the CM 
underscore the sensitivity of the B-parameter in 
(2) to changing surface roughness and also show 
that: 1) the B-parameter is also sensitive to wind 
speed, 2) under some conditions the exponent n 
differs from a value of 1.0, and 3) the value of B 
and n depend on the choice of reference levels at 
which the ambient temperature T a and wind speed 
are measured. These simulations suggest that val- 
ues of B and n can be generated for a wide range 
of conditions using the CM. 

Under minimal advection conditions, T 0 and 
T a are closely related. Above the planetary surface 
layer, however, T a and wind speed are largely 
independent of local terrain influences, are rela- 
tively homogeneous in space and therefore inde- 
pendent of small-scale horizontal variations in T 0 . 
Therefore, we will discuss the values of B and n at 
50 m. Advantages of adopting a level above the 
surface layer in the use of (2), aside from that of 
spatial homogeneity, are as follows: First, the laws 
of similarity used in calculating surface fluxes break 
down when the measurement level lies with a 
factor of about 10 times the roughness height 
above the level of the displacement height. Thus, 
over rough surfaces, such as some vegetation 
canopies (notably trees), the 2 m temperatures 
may be inapplicable to the B-method. Second, 
there are situations in remote sensing where 
screen-level measurements are not available. Third, 
50 m elevation corresponds more closely than 
screen level to the lowest atmospheric level in 
weather forecasting and climate models. 


Method for Obtaining the Parameters B and n 
and B' and n' 

The reason for using a boundary layer model is to 
generate a range of meteorological conditions, sur- 
face temperature, wind speed, surface roughness 
length, and vegetation cover from which we can 
determine the constants in (2) and (3). Since the 
results are not highly sensitive to the choice of 
vertical temperature and dewpoint profiles above 
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the surface layer, the reference set of initial atmo- 
spheric conditions for the model was chosen as 
that used by Taconet et al. (1986) in their simula- 
tions for an agricultural (wheat-growing) region in 
central France (the Beauce) during mid-July (11 
July 1983). We regard these initial conditions as 
representative of a typical sunny day during sum- 
mer in a temperate climate. It should be pointed 
out, however, that the choice of a minimum sto- 
mata] resistance in the model appropriate to wheat 
somewhat limits the results for vegetation to simi- 
lar types of crops. 

Net radiation, surface temperature, evapo- 
transpiration, and surface sensible heat flux were 
calculated as cumulative amounts over 24-h cycles 
by varying a parameter called the moisture avail- 
ability from 0.04 to 1.0 in a cycle of seven steps. 
Since moisture availability strongly modulates the 
surface temperature, this range of moisture avail- 
ability allows us to generate a large range of sur- 
face temperature. Moisture availability, as defined 
by Carlson (1986), is simply the ratio of evapotrans- 
piration to potential evapotranspiration at the sur- 
face temperature T 0 . Moisture availability also 
constitutes the link -between the surface relative 
humidity and the substrate water content, being 
set equal to the ratio of actual substrate water 
content to that at field capacity for the soil (0.34 
by volume). 

In order to obtain sufficient data to generate 
values of the model parameters, each moisture 
availability cycle was run for five surface roughness 
lengths and wind speeds and three leaf area in- 
dices, for a total of 75 moisture cycles (Table 1). 
For each moisture value, the 1300 LST surface 
temperature, air temperature (at 2 and 50 m), and 
wind speed (at 6.4 and 50 m) and the cumulative 
evapotranspiration and net radiation were com- 
puted. The rate of morning increase in tempera- 
ture, between 0800 and 1300 LST (expressed in 
°C per hour), was also calculated. We chose to 


Table 1. List of Parameters Varied in the Model and Their 
Values 


Z 0 ( m ) 

h(m) 

d 0 (m) 

Wind Speed a 

LAI 

0.01 

0.077 

0.051 

0.25 V r 

0 

0.065 

0.50 

0.335 

0.5 V r 

2 

0.13 

1.00 

0.67 

1.0 v f 

4 

0.26 

2.00 

1.34 

15 V r 


0.52 

4.00 

2.68 

2.0 V, 



“ y.-reierence wind speed. 


reduce the number of computer simulations by 
making thermal inertia a weak function of mois- 
ture availability. This expedient is justified by our 
own experience that small changes in thermal iner- 
tia do not significandy influence the results, partic- 
ularly over vegetation. 

Table 1 summarizes the roughness, wind speed, 
and leaf area index categories used in the simula- 
tions. Note that wind speed is simply stated as a 
fraction of a base (reference) wind speed for the 
11 July winds. Thus, 0.25 V T means that all the 
wind speeds in the reference vertical wind speed 
profile (from the surface to the top of the sound- 
ing) were multiplied by 0.25. 

Vegetation height ( h ) and displacement height 
( d 0 ) were made a function of surface roughness 
according to the commonly accepted relationships 
of Monteith (1975), 

z 0 = 0.13 h, d 0 = 0.67 h. 

For bare soil, the surface roughness length was also 
varied over the same range of values, although we 
recognize that the larger values of roughness would 
not likely apply to a natural bare soil cover. Air 
temperatures at 2 m, however, were not calculated 
for the 0.52 m roughness category, since the im- 
plied vegetation height would be above 2 m. 

RESULTS 

Figure 1 shows the results for seven different 
values of moisture availability plotted on log-log 
coordinates. The values of n and B are, respec- 
tively, the slope of the line and the ordinate inter- 
cept. The “best fit” straight line is defined as that 
which minimizes the residuals in the linear root- 
mean-square difference between R n2 4 — LE U (the 
surface sensible heat flux H u ) calculated from the 
model (the points in Fig. 1) and that determined 
from the straight line fit through successive pairs of 
points. 

We found that the results varied insignificantly 
with changing leaf area index in the range 2-4. 
Consequently, in subsequent discussion, we make 
no distinction between different leaf area indices, 
preferring to treat only two different types of 
surfaces: bare soil and vegetation. 

In general, simulated data adhered very closely 
to the straight line model in log-log coordinates. A 
characteristic of the best fit lines (such as that in 
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<T 0 -T a ) ‘C 


Figure 1. Values of R„ 2i - L e vs surface air temperature 
differences (T 0 - T a ) generated from a boundary layer model 
for seven values of moisture availability. The three sets of 
values, for which the best-fit straight lines are shown, pertain 
to a surface roughness (j„) of 0.065 m and a base (reference) 
wind speed for each of three surface categories: leaf area 
indices of 2.0 (A) and 4.0 (O) and bare soil (O). 

Fig. I), however, is a tendency for the lowest two 
or three points to deviate systematically from a 
straight line, although not always in the same 
sense; usually, however, the deviation was such 
that the generated values of R n24 - ££<24 fell be- 
low the line. The simulated data points using (3) fit 
a straight line with slighdy less fidelity than for (2), 
but with similar deviations at the lower end. To 
some extent, this deviation is an artifact of the 
log-log plot, which exaggerates error on the lower 
part of the scale. However, while not constituting 
a serious error, this deviation is felt to be real and 
will be discussed later in this paper. Nieuwenhuis 
et al. (1985; Fig. 3) also show a tendency for the 
data to deviate markedly from a straight line fit (in 
linear coordinates) for large values of evapo- 
transpiration (low values of R n2r 4 — LE 2 4)- 

B and n for Wind Speed and Air Temperature 
Measured at 50 m 

Results are presented in the form of analyses of B 
and n as a function of log z 0 and the 50 m wind 
speed for air temperatures measured at 50 m (Figs. 
2 and 3). Values of B [Figs. 2a) and 3a)] are close 
to 0.025. This is similar to values proposed by 



Figure 2a). Distribution of the constant B in Eq. (1) gener- 
ated by the boundary layer model for the bare soil case. The 
graph shows the constant as a function of surface roughness 
(z 0 ) and wind speed for an air temperature measured at 50 m 
and a wind speed at 50 m. 



Figure 2b). Distribution of the constant n in Eq. (1) as in 
Fig. 2a). 
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Figure 7a). Same as Fig. 2a) but for vegetation and an air Figure 7b). Same as Fig. 2b) but for vegetation and an air 
temperature at 2 m; wind speed at 6.4 m. temperature at 2 m, wind speed at 6.4 m. 


Seguin and Itier (1983). Although there is some 
nonlinearity in the distribution of B, this coeffi- 
cient generally decreases with increasing rough- 
ness and wind speed over bare soil [Fig. 2a)] and 
increases with increasing roughness and wind 
speed over vegetation [Fig. 3a)]. Differences be- 
tween bare soil and vegetation are due to the 
differing methods of calculating sensible heat flux 
from leaves and bare soil. Within the vegetation 
layer, sensible heat flux is obtained as a residual 
between the transpiration and the net radiation 
absorbed by the vegetation, whereas the bare soil 
system of equations involves a simultaneous solu- 
tion between sensible and latent heat fluxes. 

The trend in the exponent n with changing 
wind speed and roughness is similar for bare soil 
and vegetation [Figs. 2b) and 2c)]. For all practical 
purposes, however, the value of n is equal to 1.0, 
which is commonly accepted for the B-method. 

B' and n' for Wind Speed Measured at 6.4 m 
Figures 4a) and 5a), showing B', are analogous to 
Figs. 2a) and 3a), but for the temperature ten- 


dency variation of the B-method (3) with a wind 
speed at 6.4 m. The trends in B' with changing 
surface roughness and wind speed are similar to 
those in Figs. 2a) and 3a), except that the values 
differ because the units are different. The expo- 
nent n' varies rapidly for both bare soil and vege- 
tation as a function of surface roughness and wind 
speed [approximately from 1.0 to 2.5; Figs. 4b) 
and 5b)]. 

B and n for Wind Speed Measured at 6.4 m and 
Air Temperature Measured at 2 m 
Figures 6 and 7 are analogous to Figs. 2 and 3 but 
for the 6.4 m winds and 2 m air temperatures. The 
most significant differences in B between the 50 m 
and low-level values is that the former are some- 
what larger than the latter [compare Figs. 2a) and 
6a)]. However, the exponent n decreases from the 
lower to higher levels [compare Figs. 2b) and 6b)], 
thereby compensating for the increase in B with 
height. Because of this compensation, there is little 
difference in the computed values of R n24 — LE 2i 
between the 2 and 50 m results for bare soil. 
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For vegetation however, the differences be- 
tween the 50 m and lower-level temperature and 
wind levels are significant. The parameter B is 
about twice as large as the 50 m values for vegeta- 
tion [0.04-0.11 versus 0.03; compare Fig. 3a) with 
Fig. 7a)], although the variation with surface 
roughness and wind speed is similar at the two 
levels. The exponent n, while very close to 1.0 for 
both levels [Figs. 3b) and 7b)], is slightly less than 
1.0 for the 50 m results. These 2 m temperature 
values of B are similar to those proposed by Jackson 
et al. (1977). 

COMPARISON WITH FIELD 
MEASUREMENTS 

We have made an attempt to compare our results 
with local micrometeorological measurements 
made at an agricultural site operated by the De- 
partment of Meteorology at Pennsylvania State 
University. The measurements, which were made 
over wheat both, before and after harvesting in 
July, were obtained as part of a field program, 
which was in its first full summer of operations 
during 1986. The site was located in a 15 km wide 
valley about 1 km north of the foot of a mountain 
range, whose highest elevations extend about 
300-400 m above the valley floor. Air can move in 
a long fetch from the west parallel to the moun- 
tains and reach the sensors in the instrument array 
after a trajectory of about 200 m over the wheat. 
For wind directions other than westerly, the fetch 
is not as uniform, and the air moves a shorter 
distance over the wheat and may traverse other 
types of row crops. 

Comparisons were made between results of 
the measured and simulated integrated sensible 
heat fluxes. Differences between model and mea- 
surement were systematic, although measured val- 
ues of B were about one-third to one-fifth those 
obtained with the model. We feel that this discrep- 
ancy arises from errors in the measured surface 
sensible heat fluxes, which were too small. 

It is easy to demonstrate the deleterious effect 
of small systematic errors in the individual surface 
sensible heat fluxes on the 24-h averages. For 
example, consider an innocuous (but systematic) 
error of 20 W m “ 2 in surface sensible heat flux on 
a day when the maximum heat flux is 100 Win' 2 . 
This error represents an uncertainty of more than 


a factor of two in the 24-h integrated surface 
sensible heat flux. Such susceptibility of the inte- 
grated sensible heat fluxes to small error may be 
the reason why no one has yet published a statisti- 
cal comparison between measurements and mod- 
eled fluxes obtained from the B-method using con- 
stants derived from a boundary layer model. 

SOURCES OF ERROR IN MODELED 
ESTIMATES OF B, B\ n, AND n' 

Seguin et al. (1985) refer to a ±0.1 cm day -1 (or 
about 20%) error in evapotranspiration as being 
acceptable for the B-method, given the present 
state of the art. In what follows, the magnitude of 
the error pertains to the quantity R n24 — LE U and 
an error of less than 0.1 cm is considered to be 
unimportant. Errors in using (2) and (3) from this 
paper can arise from two important sources: those 
due to the model and those due to fitting the 
model results in a straight line. 

Errors in the boundary layer model are difficult 
to assess. Complex models such as the CM contain 
many different types of parametrizations. Al- 
though this model has been tested under a variety 
of conditions, it is probably difficult to verify the 
fluxes to within an uncertainty of about 20% of the 
maximum value for the day. As we pointed out, 
this magnitude of error is highly destructive of 
daily average fluxes. 

Errors in determining B and n can occur as 
the result of a deviation between model and 
straight line. With regard to Figure 1, the gener- 
ated data points diverge from a straight line at 
small surface-air temperature differences. Our re- 
sults show that virtually all the residuals in this 
region of the graphs were less than 0.1 cm. Excep- 
tions were at low surface-air temperature differ- 
ences (less than 4°C) and high wind speeds (greater 
than 4 ms' 1 ), in which errors of up to 50% were 
found. These errors become larger than 0.1 cm for 
lower wind speeds and higher surface-air temper- 
ature differences as the surface roughness in- 
creases. For small roughness lengths (less than 0.1 
m), the error in fitting a straight line was less than 
0.1 cm over the entire range of wind speeds and 
surface-air temperature differences that were cal- 
culated. 

Deviations from the best fit straight line are 
analyzed in Figure 8, which shows that the lowest 
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Figure 8. LAI * 4.0; z 0 * 0.13; base wind. * — • 2 m mid- 
summer; X — X 50 m midsummer; ® — « 2 m late summer; 
X — X 50 m later summer; . . . best fit line from 3 to lowest 
(To-T a ). 

three points suggest a straight line with a steeper 
slope and a lower intercept than the best fit fine. 
When only the net radiation in the model was 
decreased, commensurate with a change in ini tial 
conditions to late summer (September), the best-fit 
lines resemble those for the lowest three points in 
the July case. However, the lowest three late- 
summer points also deviate in the same sense as 
for the July conditions, suggesting the radiative 
forcing near the surface is an important factor in 
the determination of B and n. The same trend was 
noted for B' and n\ We believe that the weaker 
radiative forcing and smaller surface-air tempera- 
ture differences involve weaker static stability near 
the surface and therefore a change in the calcu- 
lated surface sensible and latent heat fluxes as 
dictated by similarity theory. 

In order to assess the potential error from 
choosing an incorrect B or n, we estimated a 
maximum possible error that would result from 
choosing a fixed value of B and n (or B' and n') 
without regard to wind speed and surface rough- 
ness. Let us define this maximum error as the 
maximum variation in R n24 - LE U over the entire 
domain of Figs. 2-7. This domain varies from 1 to 
9 ms -1 in wind speed and from 0.01 to 0.52 m in 
surface roughness. Using the values of B and n (or 
B' and o') at these extrema (the four comer points 


Table 2. Representative Values for B, n, B\ and n' and 
Their Maximum Variations in the Domain of z 0 and 
Winds peed 



Bare Soil 

Vegetation 


B 

n 

B 

n 

Windspeed at 50 m, 





T a at 50 m 

0.023 

1.0 

0.032 

1.0 

Max variation (cm) 

0.06 


0.37 


Windspeed at 6.4 m. 





T a at 2 m 

0.032 

1.1 

0.032 

1.0 

Max variation (cm) 

0.30 


0.63 



B 1 

n' 

B' 

n' 

Windspeed at 6.4 m. 





IT/ln 

0.062 

1.6 

0.09 

2.1 

Max variation (cm) 

0.32 


1.03 



of Figs. 2-7), we calculated the maximum range in 
R„ 24 ~ LE 24 over the domain of these figures. This 
maximum range is listed in Table 2 in the row 
labeled “max variation.” 

If one were to choose B and n (or B' and n') 
at the midpoints of these figures (where the sur- 
face roughness is about 0.1 m and the wind speed 
about 5 ms -1 ), instead of at the extrema, the 
maximum deviation of B and n from this central 
point would be about one-half that for the maxi- 
mum range across the entire domain of the figures. 
The average deviation of the integrated evapo- 
transpiration from those in which B and n are 
chosen in the middle range would be about one- 
quarter of the maximum possible error listed in 
Table 2. Dividing the values listed in Table 2 by 4 
therefore gives us an estimate of the average error 
in using (2) and (3), if B and n (or B' and n') 
were chosen from the center of these graphs. This 
error is generally less than 0.1 cm in most, but not 
all, circumstances. 


CONCLUSIONS 

The present study addresses the variation in the 
parameters B and n in the B-method. These pa- 
rameters are sensitive to wind speed and rough- 
ness, especially over vegetation. Values of B are 
similar to those of Jackson et al. (1977) for vegeta- 
tion and Seguin and Itier (1983) for bare soil. It is 
interesting to note that Rambal et al. (1985) also 
find that B increases with increasing vegetation 
height and density. For n = 1.5 Rambal et al. 
obtained B = 0.032 for shrubland and B = 0.068 
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for middense woodland, values which resemble 
ours for bare soil and vegetation, respectively. For 
all practical purposes, however, we find that the 
exponent n was 1.0 for both the bare soil and 
vegetation cases. 

We suggest that the B-method might be more 
generally applicable to regional-scale remote sens- 
ing if the reference level for air temperature and 
wind speed were taken at 50 m, rather than at 2 m 
because the temperature and wind speed at the 
50 m level are relatively insensitive to local surface 
inhomogeneities. For this reason we show the dis- 
tribution of B and n for air temperatures mea- 
sured both at the customary 2 m elevation (6.4 m 
for wind speed) and at 50 m. In view of all the 
errors inherent in the B-method, it may be suffi- 
ciently accurate to use regional-scale surface wind 
and air temperature measurements and to stratify 
the wind speed, surface roughness, and vegetation 
type into just two or three categories. 

We also propose using the morning rise in 
radiometric surface temperature [as suggested by 
Wetzel et al. (1984)] in the B-method, because this 
does not require knowledge of an air temperature. 
Although the data for the morning rise in surface 
temperature fit straight lines in log-log coordi- 
nates, the coefficient B' and the exponent n' are 
highly sensitive to wind speed and surface rough- 
ness. Moreover, the latter differed considerably 
from unity for most combinations of wind speed 
and surface roughness. 
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ABSTRACT 

Lynn, B.H. and Carlson, T.N., 1 990. A stomatal resistance model illustrating plant vs. external control 

of transpiration. Agric For. Meteorol. , 52: 5-43. 

A stomatal resistance model is proposed that is specifically designed for use within atmospheric 
boundary layer models. Stomatal resistance is expressed as a product of linear discontinuous func- 
tions of leaf water potential and solar flux; vapor pressure deficit is included indirectly by allowing 
for a gradient of leaf water potential between the surface and the interior of the leaf. An important 
assumption is that stomatal resistance increases sharply beyond a critical value ofleaf water potential 
and/or critical solar flux. The model is able to simulate a transpiration plateau; this plateau may occur 
because of high atmospheric demand and/or a limitation in soil water. 


INTRODUCTION 

Evapotranspiration is highly dependent upon the resistance exerted by the 
plant on the flow of water from soil to atmosphere. Owing to the complexity 
of calculating plant resistance, one approach has been to represent the vege- 
tation canopy as a uniform “big leaf’. This method employs a combination 
equation, such as the Penman-Monteith formula (Monteith, 1975 ), in which 
the soil/plant (canopy) resistance constitutes a basic unknown. More com- 
plex models allow for the calculation of both a soil and vegetation flux, which 
depend, respectively, on soil and plant (stomatal) resistances (DeardorfT, 
1978; Soer, 1980; Dickinson, 1984; Halldin et al., 1984/85; Lindroth and 
Halldin, 1986;Taconet etal., 1986a, b; Sellers and Dorman, 1987; Wetzel and 
Chang, 1987;Lhomme, 1988; Abramopoulos et al., 1988). 

There is a general consensus among plant scientists that stomatal resistance 
is correlated with a number of factors. These are short-term changes in leaf 
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water potential, significant drying of the soil (very negative soil water poten- 
tial), vapor pressure deficit, solar flux, leaf temperature, and ambient carbon 
dioxide (Turner, 1970, 1974; Van Bavel, 1974; Jarvis, 1976; Tan and Black, 
1976; Takami and Uchijima, 1977; DeardorfT, 1978; Farquhar, 1978; Fed- 
erer, 1979; Takami and Yukimura, 1979; Singh and Szeicz, 1980; Fisher et 
al., 1981; Farquhar and Sharkey, 1982; Kaufmann, 1982; Choudhury, 1983; 
Zeiger, 1983; Dwyer and Stewart, 1984; Halldin et al., 1984/85; Avissar et 
al., 1985; Choudhury and Idso, 1985;Komer, 1985; Lindroth, 1985; Simpson 
et al., 1 985; Grantz and Zeiger, 1 986; Lindroth and Halldin, 1 986; Baldocchi 
et al., 1987; B. Acock, personal communication, 1989). It has also been pro- 
posed that there is an additional dependence of stomata! resistance on soil 
water potential, caused by the transport of cytokinins from the root in re- 
sponse to soil drying (Gollan et al., 1986; Schulze, 1986; Munns and King, 
1988). 

In order to express the relationship between these factors and stomatal re- 
sistance, it is necessary to combine them in a deterministic manner. There are 
numerous examples in the literature in which stomatal resistance is expressed 
as either an additive or multiplicative model. In some cases, deterministic 
expressions are both additive and multiplicative (Table 1 ). Almost all models 
are expressed as a function of solar flux, soil or bulk mesophyllic leaf water 
potential, and vapor pressure deficit (Table 2). 

The purpose of this paper is to formulate a stomatal resistance model that 
also incorporates an analytical solution for leaf water potential. A unique as- 
pect of this model is that it allows for epidermal control of leaf water potential 

TABLE 1 


Stomatal resistance functions by type 


Model 

Type I 

Type II 

Type III 

Avissar et al. (1985) 


+ 


Baldocchi et al. (1987) 


+ 


Choudhury (1983) 


4- 


Choudhury and Idso ( 1 985 ) 


+ 


DeardorfT (1978) 



+ 

Dwyer and Stewart ( 1984) 

+ 



Federer (1979) 

+ 



Halldin etal. (1984/85) 


+ 


Jarvis (1976) 


+ 


Kaufmann (1982) 



+ 

Simpson et al. (1985) 



+ 

Singh and Sziecz (1980) 



+ 

Takami and Yukimura (1979) 


+ 



Each equation has been assigned a type, i.e. additive (Type I), multiplicative (Type II), while 
Type III is a hybrid of Type I and Type II. 
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TABLE 2 


Siomatal resistance functions by variable 


Modeler 

/(wat) 

AS) 

A n 

AT) 

Avissaret al. ( 1985) 


s 

I’a 

r, 

Baldocchiet al. (1987) 

V\ 

PAR 

v\ 


Choudhury ( 1 983 ) 



* 

* 

Choudhury and Idso ( 1985) 

W\ 

5 

* 

* 

DeardorfT ( 1978) 

0vo/<9 v 

5 

5 

* 

Dwyer and Stewart (1984) 

0v 

5 


* 

Federer (1979) 

Wi 

S 

v. 

T a 

Halldin et al. (1984/85) 

- 

s 

K 

* 

Jarvis ( 1976) 

¥\ 

s 

v. 

T t 

Kaufmann ( 1982 ) 

* 

PAR 

V, 

* 

Simpson et al. (1985) 

y'xb 

* 

K 

* 

Singh and Sziecz ( 1976) 

* 


* 

* 

Takami and Yukumura ( 1979) 

V\ 

* 

* 

* 


Symbols: (*) indicates not included and (-) is not applicable. Table 2 is a listing of model types 
and the parameters included in each./(wat) refers to parameters relating r, to water content of 
the soil or leaf. The other headings pertain to the mathematical functions governing r, for solar 
flux (/(S) ), vapor pressure deficit (/( V ) ) and temperature (AT)). R ncl is the net radiation, 5 
is the incident solar flux, and PAR is the photosynthetically active radiation (proportional to 
0.47 S (van Bavel, 1974)). % is the leaf water potential, is the pre-dawn (base) xylem 
potential, and 0 VO and 0 V refer to the surface and root zone soil water content. Land Tare the 
vapor pressure deficit and temperature, respectively. Subscript a refers to the atmosphere while 
1 refers to the leaf or leaf-air difference. 


by specifying a gradient of leaf water potential between the surface and the 
interior of the leaf. This gradient is dependent on the vapor pressure deficit. 
Leaf water potential is calculated as a function of soil water potential, solar 
flux, vapor pressure deficit, leaf boundary layer resistance and root-stem re- 
sistance, thus incorporating indirectly the effect of these parameters on sto- 
matal resistance. 

An assumption in the model is that there exists a critical value of leaf water 
potential at which stomatal resistance increases sharply with decreasing leaf 
water potential. A number of researchers have justified this assumption 
(Turner, 1 974; Boyer, 1 976; Thomas et al., 1 976; Denmead and Millar, 1976a; 
Kateiji, 1979; Pospisilova and Solarova, 1980; Idso, 1983; Dwyer and Stew- 
art, 1984; Schulze et al., 1987). However, in this paper we assume that the 
critical water potential refers to that of the epidermis as suggested by Sheriff 
(1984) and Schulze (1986). 

In order to allow for the interaction between soil, plant, and atmosphere, 
we have incorporated the stomatal resistance model within a canopy model 
(which is an extension of the model by Taconet et al. ( 1 986a) (Fig. 1 ) ), and 
a one-dimensional boundary layer model of Carlson (1986) and Taconet et 
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Fig. 1 . Structure of the plant-canopy model, showing the exchange of sensible and latent heat 
fluxes between plant, atmosphere and surface (see text for details). 

al. ( 1 986b ) . In this paper, we discuss the stomatal resistance formulation and 
the solution to the leaf water potential equation, and discuss its significance. 
We then show how the model is consistent with some field measurements. 

PLANT MODEL 

The flow of water from soil to the atmosphere can be described by a flow 
diagram as seen in Fig. 1 . Table 3 contains a listing of all the variables in this 
figure and the equations below. 

Water flow from soil to leaf 

The movement of water from the soil to the leaf begins in the soil root zone 
with soil water content 6 V or soil water potential y/ v Following Ohm’s law the 
water flux from soil to leaf is given by 

F v , = {y/ i -y/ ] -H)/{Z p +Z t ) (1) 

where: the variables are the soil water potential (y/ s , bar), and the mesophyl- 
lic leaf potential (y/ h bar); H is the gravitational potential (bar); Z p is the 
root-xylem resistance (bar (W m~ 2 ) ~ 1 ); the soil root interface resistance is 
Z g (bar (W m -2 ) -1 ). Here H=p„gh, where:/)* is the density of liquid water 
(kg m -3 ); g is the gravitational constant (m s -2 ); h is the average height of 
the leaves above the root zone (m). Equation ( 1 ) has gained wide acceptance 


STOMATAL RESISTANCE MODEL 


9 


with plant scientists (Cowan, 1965 (who cites a seminal paper by van den 
Honert); Boyer, 1971;FeddesandRijtema, 1972; Camacho etal., 1974; Den- 
mead and Millar, 1976b; So et al„ 1976; Meyer and Green, 1980; Katerjl and 
Hallaire. 1984; Passioura, 1984;PassiouraandMunns, 1984; Choudhury and 
Idso, 1985). 

Description of soil water variables 

In the boundary layer model, soil water content is initialized with the mois- 
ture availability parameters of the surface (M c ) and of the root zone (M). 
M 0 and M are equated to their respective soil water contents by the somewhat 
arbitrary but empirically justifiable assumption that M 0 = 9 vo /6 vc and M= 0 V / 
0 VC . 0 VO and 0 V are the soil water contents of the surface layer and root zone, 
respectively, and 0 VC is the value of soil water content at field capacity. 

As noted in eqn. ( 1 ), the flow of water from soil to leaf is dependent upon 
soil water potential. Soil water potential (bar), for a given soil type, is 
related to 0 V (cm 3 cm -3 ) using a parametric equation discussed in Campbell 
(1974). This equation is: 

V'g = VA«(0v/0vsal) _& 

where: is the value of soil water potential at saturation (bar); 0 v »t is the 

value of soil water content at saturation (cm 3 cm -3 ); b is a constant. 

Z g expresses the resistance of the flow of liquid water from the soil to the 
roots, in response to a gradient of water potential that increases as the soil 
dries (Cowan, 1965). Various ways exist to determine Z g , e.g. Cowan (1965), 
Feddes and Rijtema (1972), Denmead and Millar (1976b) and Federer 
(1979, 1982). The choice for the present model is based on a scheme of 
Choudhury and Idso (1985) which was originally proposed by Feddes and 
Rijtema (1972). They let 

Z g = (0.0013/(Z efTJ K g ))/: 1 (3) 

where: 0.001 3 (m 2 ) is the ratio of a parameter relating root distance and ge- 
ometry to the reciprocal of the effective rooting depth; Z eff is the effective 
rooting depth ( m ) ; K s is the soil hydraulic conductivity ( m s - 1 ) , and is gi ven 
by Campbell (1974) as: 

K t = KM/e v „ x ) 2b+2 (4a) 

or, 

r» M /v' t ) 2+2/b < 4b > 

where: is the maximum conductivity at field saturation of the soil (m 

s -1 ); 0 vsat is the saturation volume water content; b is a constant that is de- 
pendent on soil type; y/ gsat is the value of the hydraulic potential at field satu- 
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TABLE 3 


List of variables and coefficients 


Variable Comments 

a leading coefficient, eqn. ( 7 ) , set equal to ^min 

b coefficient in eqns. (2), (3), (4) 

b the slope of the subcritical part of thc/(^ e ) function; 


Units 
s m* 1 

s m~ 1 bar -1 


b 2 

P 


Cl 

C2 

Cp 


<?af 

e x (T x ) 

F w 

f(C0 2 ) 

/(vO 

AS) 

ATi) 

g 

y 

h 

H 

H t 

K x 

U 

L*E f 

LPSI 

M 0 

M 

P 

Pw 

r. 


r*> 

'd 
r cut 

r\ 

^min 


W t >W c 

the slope of the supercritical part of the/(^ c ) function; 

constant describing the difference between mesophyllic and 
leaf epidermal water potential which is divided by vapor 
pressure deficit 

the slope of/(S) between S 0 and S where S> S c 
the slope of AS) between S c and S where 5<5 C 
the specific heat of dry air at constant pressure 
& ¥ = 0 (y/ e > 1 and ip c =y/ c in the second 

term on the righthand side of eqn. (8a) 

<5 S =0 (S>S C )\ <5,= 1 {S^So), and S=S C in the second term 
on the righthand side of eqn. ( 8b ) 
leaf-air boundary vapor pressure 

the saturation vapor pressure at the temperature of the leaf 
flow of water from soil to leaf 
stomatal resistance factor due to carbon dioxide 
stomatal resistance factor due to epidermal leaf water 
potential 

stomatal resistance factor due to solar flux 
stomatal resistance factor due to leaf temperature 
gravitational constant 
/>C p /0.622I e 

the average height of the leaves above the root zone 

gravitational potential 

foliage sensible heat flux 

the soil hydraulic conductivity 

a conversion factor 

latent heat of vaporization of water 

foliage latent heat flux 

Leaf Potential Stress Index 

A/ o =0 V q/0 vc (surface moisture availability) 

M=ejQ y c (moisture availability of the root zone) 
the atmospheric pressure 
the density of liquid water 
air resistance in surface layer 

the resistance for heat and water vapor flux for the interleaf 

air spaces . 

air resistance between the ground and the interleaf air 

spaces 

air resistance in transition surface layer 
“critical” resistance 
cuticular resistance 
leaf resistance 

minimum stomatal resistance 


s m - ' bar -1 
bar mbar” 1 


sm- 1 (Wm- 2 )- 1 
s m“ 1 (W m“ 2 ) 1 
Jkg-'K ’ 1 


mbar 

mbar 

Wm' 2 


m s' 2 
mbarK“ l 
m 
bar 

Wm- 2 
m s~* 

bar (Wm- 2 )*' 

Jkg- 1 

Wm" 2 

bar 


mbar 
kgm“ 3 
s m” 1 

s m _l 

s m” 1 
s m _1 
s m* 1 
s m” 1 
s m" 1 
s m~ l 



STOMATAL RESISTANCE MODEL 


11 


TABLE 3 (continued) 


Variable 

Comments 

* 

Units 

r. 

stomatal resistance 

s m~ 1 

6 K 

soil water content of the root zone 

cm 3 cm" 3 

0VC 

soil water field capacity 

cm 3 cm" 3 

0v„ 

surface soil water content 

cm 3 cm -3 

^vsat 

value of soil water content at saturation 

cm 3 cm" 3 

Vc 

critical leaf water potential, which is the intersection of the 
two lines with slopes b\ and b 2 

bar 

Vc 

epidermal leaf water potential 

bar 


soil water potential 

bar 

v* 

critical soil water potential 

bar 

Vi 

mesophyllic leaf water potential 

bar 


the value of soil water potential at saturation 

bar 

s « 

the solar radiation threshold, which is at the intersection of 



lines formed by the slopes c \ and c 2 

Wm" 2 

So 

a maximum solar radiation value, loosely defined as the 



value at which light saturation of the leaves is reached 

Wm' 2 

(7 

pL t 

kgm" 3 Jkg" 1 

T. 

air temperature of the surface layer 

K 

T, t 

temperature of the interfoliage air spaces 

K 

T t 

temperature of the ground surface 

K 

r, 

temperature of the leaf surface 

K 

V 

the difference between e t { T x ) and e, f 

mbar 

WPSI 

Water Potential Stress Index 

bar 

Z'fT 

the effective rooting depth 

m 

Z t 

the soil root interface resistance 

bar (Wm" 2 )" 1 

z p 

root-xylem resistance 

bar (Wm- 2 )- 1 

Z f00 | 

root resistance 

bar (W m -2 )- 1 

7 

stem 

stem resistance 

bar (W m“ 2 )- 1 

z, 

Z g -1- Z p 

bar (Wm” 2 )* 1 


ration. K\ is a conversion factor which enables Z g to be expressed in bar (W 
m -2 ) -1 , given as 0.4x 10 -10 . 

y/ t and Z g (AT g ) vary with soil moisture (6 V ) in a parametric manner as 
determined by 0 vsat , 1 //^,, K^, and b which are all functions of soil type. 
becomes very negative, and Z g large (AT g small) when 6 V approaches a critical 
value. This critical value may sometimes be defined as the “nominal wilting 
point”, or the soil water content at which the leaves have lost or lose turgor in 
the early morning. The values of these parameters can be found in Cosby et 
al. ( 1984). Those for loam were found to agree with measurements made at 
a local agricultural field site (Rock Springs), provided that 0 VS at> which is 0.43 
as proposed by Cosby et al. (1984), is replaced by 0.34, which is the field 
capacity at Rock Springs. Henceforth, this value of field capacity will be used 
as a basis for calculation throughout this paper. 
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Root and stem resistance 

The internal plant resistance, Z p , consists of the root and stem resistance, 
referred to as Z rooi and Z stem , respectively (see Fig. 1 ). In principle, Z p can be 
determined by dividing the difference between a value of leaf water potential 
and the pre-dawn leaf water potential by the transpiration rate at the time at 
which the leaf water potential is made. This calculation is often complicated 
by a hysteresis effect produced by water storage in the plant (Kateiji and 
Hallaire, 1984). 

Unfortunately, very few measurements of Z p have been published, proba- 
bly because of the difficulty in measuring leaf water potential and transpira- 
tion simultaneously. Some estimates of Z p can be obtained from, or are listed 
in Boyer (1971), Denmead and Millar (1976b), Jarvis (1976), Running 
(1980a); Kateiji et al. (1983), and Abdul-Jabbar et al. (1984). Units may 
be expressed in s, h, or in bar sera' 1 ; but here, the preferred units are in terms 
of the drop in potential from root to xylem divided by the flux per unit leaf 
area ( bar ( W m * 2 ) ~ 1 ) . The conversion from s to bar ( W m “ 2 ) “ 1 is given by 
the factor l/(/? w L e x 10), which equals 0.4x 10 -10 . (The value of 10 is the 
conversion from m to bar (10m bar" 1 ) . ) The conversion from bar s cm ~ 1 
to bar (Wm~ 2 ) -1 is given by the factor 1 00/ (/x*X e ), which equals 0.4 x 10 -7 . 
A typical value of Z p (obtained from Kateiji and Hallaire ( 1984) ) is 0.047 
bar (Wm- 2 )- 1 . 

It should be noted, however, that the value of Z p varies by species (Abdul- 
Jabbar et al., 1984). Z p also varies during the growing season because it is a 
function of the number, width, age of the xylem, and the root mass (B. Acock, 
personal communication, 1 989 ). We suggest that Z p may vary inversely with 
crop height as is indicated by the results presented later in the paper. Z p is an 
important parameter, because sensitivity tests presented in this paper show 
that the transpiration rate and the leaf stomatal resistance are indirectly re- 
lated to the choice of Z p . 

Hailey et al. (1973), Bunce (1978), Jones (1978), Kateiji and Hallaire 
(1984), and Passioura and Munns (1984) caution that the application of 
Ohm’s laws to plants is not always applicable. For instance, Bunce (1978) 
found that transpiration continued to increase in some species, (e.g. soybeans 
and cotton), while leaf water potential remained constant. Various theories 
have been proposed to explain this phenomenon. Bunce suggests that a de- 
crease in root resistance due to root elongation may decrease the value of Z p . 
However, Kateiji and Hallaire (1984) suggest that large demand causes the 
plant to release internal sources of moisture for transpiration, an explanation 
that we favor. This would suggest that the value of Z p remains relatively un- 
changed with transpiration. Thus we assume that root-stem resistance re- 
mains constant over a range of leaf water potential. 
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Water flow from the leaf to the atmospheric boundary layer 

The flux of water vapor from the leaf to the boundary layer is given as 

L'E { = (pC p /y)V/ (f| + r af ) ( 5 ) 

where: V is the difference between the saturation vapor pressure at the tem- 
perature of the leaf (e 5 (r,)) and the leaf-air boundary vapor pressure (e af ) 
(mbar); y is equal to PC p /0.622L c ( P is the atmospheric pressure (mbar), 
C p is the specific heat of dry air at constant pressure ( J kg" 1 K~ 1 ), and L e is 
the latent heat of vaporization of water (J kg 1 ); the leaf resistance (h) (s 
m -' ) is related to the stomatal resistance (r s ) (s m -1 ) by the formula for 
parallel resistances 

l/r, = l/r s + l/r cut (6) 

where: r s is the stomatal resistance; r cu , is the cuticular resistance (s m ), r af 
is the resistance for heat and water vapor for the interleaf air spaces (s m ) 
and is a function of the friction velocity and leaf area index (LAI) (Thom, 
1972; Taconet et al„ 1986a). The latent flux is linked to the temperature and 
specific humidity of the boundary layer by atmospheric resistances r 3 and r b 
(sr 1 ). 

STOMATAL RESISTANCE MODEL 
Linear discontinuous model 

We adopt a multiplicative model which is based on the product of functions 
f(y/ e ),f(S),f{Ti), and /(C0 2 ). These are the inverse fractional conductance 
due' to leaf water potential, solar flux, temperature, and carbon dioxide, re- 
spectively. (The effect of vapor pressure deficit ( V) is treated separately, and 
is discussed below; the transport of hormones from the root is not explicitly 
accounted for. ) 

Given the assumptions made above, the equation for r s is written as: 

r s = a/( 5 )/( y/ c )/( T\ )/( C0 2 ) ( 7 > 

where a has the dimensions of resistance and is equated with the minimum 
stomatal resistance r min (defined as the resistance measured when y/ e is zero, 
the leaves are saturated by sunlight, the leaf temperature and carbon dioxide 
concentration are optimum and the vapor pressure deficit is equal to zero). 
(The reader need be aware that the definition of minimum stomatal resis- 
tance is not universally expressed in this way because a minimum measured 
value will vary with vapor pressure deficit, solar radiation, leaf (or soil) water 
potential, and is a function of plant type and plant phenology. ) All of the 
functions may vary from 1 .0 to infinity. 
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Because of the limited availability of data, a simplified representation of 
the governing factors in eqn. (7) is essential. Accordingly, we propose a “dis- 
continuous linear model” similar to that discussed by Fisher et al. (*1981 ). 
Unlike Fisher’s model, however, the exponential behavior of f{S), f(y/ c ), 
and /(C 0 2 ) is represented by a pair of straight lines (Fig. 2) whose 
intersection defines a critical value, S c , y/ Q , T c and C0 2 . Figure 2 suggests that 
an appropriate straight line fit to exponential functions can capture the essen- 
tial discontinuous nature of the functions without great loss of accuracy. 

The following is a mathematical description of the factors f(y/ e ) and f{S): 

A^) = (r min +b^ e + b 2 (^c-^c)3 v )/r mm (8a) 

where: r min is the minimum value of resistance; b x is the slope of the subcriti- 
cal part of the f(y/ t ) function; y/ e >y/ c ', b 2 is the slope of the supercritical part 
of the f{y t ) function; i// c is the critical leaf water potential, which is 

the intersection of the two lines with slopes b x and b 2 , 3^=0 (y/ c > yr c ); A = 1 
( Vc), and y/ t =y/ c in the second term on the right-hand side of eqn. (8a). 

AS) = (r min +c, (S a -S)+c 2 (S c -S)3 s )/r min (8b) 

where: C\ is the slope of f(S) between S 0 and S c ; c 2 is the slope of f(S) between 
5 C and 5 where 5'<5 C ; S c is the solar radiation threshold, which is the intersec- 
tions of lines formed by the slopes Ci and c 2 ; S a is a maximum solar radiation 
value, loosely defined as the value at which light saturation of the leaves is 
reached; <5 S =0 ( S>S C ), 3,= 1 (S^S C ), and S=S C in the second term on the 
right-hand side of eqn. (8b). 

It is possible to define the response of the stomata to temperature and car- 

I 



Fig. 2. Schematic illustration of /(^ e ) function (solid curve) with changing leaf water potential, 
starting from a base value y/ t . The vertical dotted line pertains to the critical leaf water potential 
y/ c . The sloping dashed lines show the linear approximation to the function by two straight lines 
intersecting at ^ c . Also shown are the relative values off(y/ c ) f(w), and /(^,) at a given y/ ty as 
indicated by the horizontal dashed lines. 
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bon dioxide as noted above. There is, however, relatively little data on which 
to base these functions, although the response of the stomata to these factors 
may be exponential (Jarvis, 1976; Morison, 1987; Baldocchi et al., 1987). 
Taking into account the findings of the noted authors, we assume that T\ is 
within an operational range of temperature, and that the ambient carbon 
dioxide concentration has no effect. Thus /( 71) = 1.0 and /(CC> 2 ) = 1.0. It 
should be apparent that the use of this linear formulation for the governing 
factors provides a simple means of representing the effects of each factor. 

If S=S 0 , the multiplicative model reduces to an additive model. If 5 is less 
than S 0 , then ( S-S 0 ) multiplies the effect of This suggests that an 
underlying assumption of this model is that values of solar flux less than S 0 
affect the value of leaf water potential. 

A survey of the literature suggests that larger values of y/ c correspond to 



Fig. 3. Stomatal resistance vs. leaf water potential (yO and solar flux ( S ) is shown in a three- 
dimensional plot for y/ c = — 16 bar, and S c =225 W m J . 
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temperate latitude crops (typically - 13 to - 16 bar), while the for- 
est trees ( such as a pine ) might be as small as — 22 bar ( Running, . 

trees and plants growing in a desert environment might have even lower crit- 
ical leaf water potentials, and t// c may also be a function of soil water typ e an 
water stress (Dwyer and Stewart, 1984; B. Acock, personal communication, 
1989; J. Norman, personal communication, 1989). The value of S c typica y 

^Hgure^^shows how ^the 'threshold values dramatically affect resistance. 
There are relatively small changes in resistance until either the leaf water po- 
tential reaches its critical value {y e ) or the solar flux becomes less than its 
critical value (S c ). Beyond the critical thresholds, the rate of change in r s 
most strongly influenced by b 2 and c 2 . It is worth noting that either solar in- 
tensity or leaf potential can limit transpiration, but owing to the small value 
of S c the leaf water potential is more important than S c in cons rami g 
transpiration flux. 


The effects of vapor pressure deficit 

One definition of vapor pressure deficit is the difference between the satu- 
ration vapor pressure of the air above the canopy and the actual vapor pres- 
sure of the air (similar to the dew point depression). This definition is used 
because it is relatively easy to measure and corresponds to a parameter in the 
Pemttau-Monteith equation (Monteith, 1 975 > Some 
such as Camacho et al. (1974), Ackerson and Kneg ( 1977), and Luxmoore 
et al ( 1981 ) prefer to work with the leaf-air humidity difference, w 
more fundamental quantity because it represents the actual vapor gradient 
experienced by the leaf. Since the mathematical relationship between these 
definitions of vapor pressure deficit is not unique (aoudhury and Mot.te.th, 
1986) it is more physically correct to use the definition of V that represents 
the vapor pressure difference between the leaf surface and that in the interleaf 

airspaces. We adopt this definition of V. trk 

Most researchers have characterized the response of stomatal resistance to 
vapor pressure deficit as either linear, exponential, or both, depending up 
species. The linear relationships have been P r °P°f d ^ ederer [ L s (’ 
Choudhury (1983), Halldin et al. (1984/85), Feldhake and Boyer (1985), 
Simpson et al. (1985) and Lindroth (1985). ^P°^ formulations have 
been proposed by Jarvis (1976), Tan and Black ( 1976), Avissar et al. ( 1985 ) 
and Baldocchi et al. (1987). Camacho et al. (1974), Ackerson and Kneg 
(1977) and Gollan et al. (1985) have found the stomatal resistance to be 
both linear and exponential, depending on the crop or range of vapor pressure 

de After reviewing the many stomatal resistance formulations summarized in 
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Table 1 , we were unable to understand the basis for a direct effect of vapor 
pressure deficit on stomatal resistance. It may be possible to explain the re- 
sponse of the stomata to F as a response of the guard cells to epidermal leaf 
water potential. At the same time, we recognize that transpiration is corre- 
lated with vapor pressure deficit and with a gradient of water potential within 
the leaf. Therefore we postulate two leaf water potentials, one for the meso- 
phyll (y/j) and one for the epidermis (^ c ) (near the guard cells) such that 

(9) 

where is a constant. Gradients of water potential between the mesophyll and 
epidermis have been measured (Shackel and Brinckmann, 1985; Shackel, 
1987; Frensch and Shulze, 1988). 

A virtue of eqn. (9) is that it captures the “feedforward” effect described 
by Farquhar (1978), Sheriff (1984) and Schulze (1986). The feedforward 
effect is associated with an increase and then a decrease in transpiration as V 
increases. This is accompanied by a continued decrease in epidermal leaf water 
potential, while mesophyllic water potential increases beyond the point of 
transpiration decrease. Komer (1985) has observed the feedforward re- 
sponse in mountain maple, silver beech, and Scots Pine. Running (1976), 
and Sellers and Dorman (1987) have also found that the stomata of other 
coniferous trees show a similar response to vapor pressure deficit. Sheriff 
( 1977), J.A. Bunce (personal communication, 1989), and Baldocchi et al. 
(1987) have noted that the response of the stomata to vapor pressure deficit 
varies among crop plants. Others who have investigated the feedforward re- 
sponse are Zieger ( 1983), Sheriff ( 1984), El-Sharkaway et al. ( 1 984), Xu et 
al. (1984), Roessler and Monson (1985), and Grantz and Zeiger (1986). 

Solution for leaf potential 

Although the representation of natural processes within the soil-plant can- 
opy model is complex, the system of equations presented above allows for a 
mathematically simple derivation for in terms of soil moisture, solar flux, 
vapor pressure deficit, interleaf air resistance, and plant internal resistance. 

Assuming that the plant capacitance or storage capacity, as defined by Jones 
( 1983), is small, then 

L t E ( =F„ 0°) 

where F w is the flux of water from the root zone to the leaves (per unit leaf 
area). Equation (10), stating that the water flux entering the roots (F w ) is 
equal to the transpiration rate assumes a “steady state (Jarvis, 1976), 

and results in a quadratic solution for y/ e , given a linear relationship between 
r s and y/ e . 

Before solving for y/ e , it is necessary to test whether y/ c is below or above the 
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threshold value in order to determine the value of 5 V . To do so, a critical soil 
water potential ( y / gc ) is defined as the minimum soil water potential that can 
meet the evaporative demand without y/ c becoming less than y/ c . This is de- 
rived from eqn. (10) by setting the leaf water potential equal to that of the 
threshold epidermal leaf water potential, at which point r s is equal to a critical 
resistance r c1 . Rearranging then yields 

V'gc = {Vc + fiV) + ( VaZ, )/{r af + [r cut r c ,/(r cut + r cl ) ]} + H (11) 

where Z t = Z p +Z g . If the critical soil water potential is less than the value of 
the soil water potential, the subcritical solution is the correct solution because 
y/ e will be greater than y/ c . If, however, the critical soil water potential is greater 
than the value of the soil potential, is less than y/ c , requiring the supercrit- 
ical solution. 

After solving for y/ gc , the solution for y/ c is obtained by solving the following 
equation, which is obtained by combining eqns. ( 1 ), (5), (6), (7), (8a, 8b), 
and (9) in ( 10) to yield: 

A^ e 2 -t-B^ e +C=0 (12a) 

The coefficients A, B, and C contain all of the independently specified or 
calculated parameters listed in these equations. This quadratic equation is 
solved by conventional algebra to yield two roots. The following is selected as 
the correct root because it is consistent with reasonable values of leaf water 
potential: 

y/ t — (B— (B 2 — 4AC) 1/2 )/2A (12b) 

(A brief derivation of (eqn. 12) is shown in the Appendix.) 

Plant stress 

The critical ground water potential y/ tc , derived in eqn. (11), describes the 
changing relationship between plant water demand and soil water supply, and 
resembles the plant water supply function defined by Cowan (1965). It is 
apparent from eqn. (11) that if V is greater than zero, then y/ gc will be greater 
than y/ c and that y/ gc increases as V increases; y/ gc also increases as y/ t decreases 
through Z g . This suggests that both V and y/ % determine the limiting value of 
soil water potential that may support transpiration without a decrease in y/ e 
below y/ c . 

It is possible to define two measures of plant stress, which quantify the de- 
pendent water relationships noted above. These are the Water Potential Stress 
Index (WPSI) and the Leaf Potential Stress Index (LPSI) given as 

WPSI = v/gc — 

and 


(13a) 
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LPSI = y/ c - y/ e ( 1 3b ) 

When WPSI and LPSI are positive, the stomata are closing at a rate equal to 
the supercritical part of the discontinuous curve. We suggest that these in- 
dices are measures of plant stress, and it may be possible to use them to cal- 
culate a cumulative effect on plant growth. 

THE BOUNDARY LAYER MODEL 

The boundary layer model is a time-dependent, initial value model that 
contains several layers for heat and moisture transport. It consists of two soil 
layers, a bare soil surface, a layer of vegetation, a surface air layer and a mix- 
ing layer. The model starts from a set of initial values and advances at ap- 
proximately 4-min time steps from near sunrise for a period of nearly 24 h. 
In coupling the plant component with the full boundary layer model, feed- 
back is allowed between the rates of evapotranspiration and sensible heat flux 
and the specific humidities and temperatures within and above the canopy. 
The reader is referred to Carlson (1986), and Taconet et al. (1986b) for 
details. 

SCALING UP FROM A LEAF TO A CANOPY 

Similarity theory suggests (Lindroth and Halldin, 1986) that the calcu- 
lated average leaf resistance (r,) needs to be divided by the LAI in order to 
scale from a leaf to a canopy. This is equivalent to multiplying the fluxes of 
sensible and latent heat by LAI. 

In this paper we consider r, (which contains both the stomatal and cuticular 
resistances) to be calculated for a single sided upfacing leaf near the sunlight 
top of the plant canopy; thus r, is not an average resistance over the depth of 
the canopy. In order to account for the extinction of S (and thus an increase 
in y/ e ) within the canopy and the interaction between plants, we also multiply 
r, (divide the sensible and latent heat fluxes) by a shelter factor to give a 
foliage flux (Taconet et al., 1 986b), which accounts for these differences in a 
highly empirical manner. The ratio LAI/P, varies gradually from a value of 
1 .0 to 2.0 for LAI varying between 2 and infinity. Thus the inclusion of a 
shelter factor restrains large variations in evapotranspiration and sensible heat 
fluxes that would occur because of large variations in LAI. 

Henceforth all fluxes calculated within the boundary layer model, and re- 
ferred to in the figures have been scaled by LAI/P and therefore refer to an 
entire plant canopy, and to a horizontal surface area. It should be noted that 
this scaling is not strictly correct because it treats r cut in the same way as r s . 



20 


B.H. LYNN AND T.N. CARLSON 


RESULTS AND DISCUSSION 

* 

This section consists of three parts. In the first and second parts, we discuss 
how moisture availability, root-stem resistance, critical leaf water potential, 
and (feedforward effect) effect various plant and atmospheric variables. To 
do so, we use the stomatal resistance model independently and within the 
boundary layer model. We also discuss an interesting phenomenon, which is 
the appearance of a plateau in the transpiration. This plateau is associated 
with a transient increase in stomatal resistance during the midday period. 
Two indices which are related to the appearance of the plateau are also shown: 
the Water Potential Stress Index (WPSI), and the Leaf Potential Stress Index 
(LPSI). The third section, pertaining to measurements made during the 
HAPEX-MOBILHY experiment, presents experimental evidence suggesting 
the importance of critical leaf water potential, and root-xylem resistance. (The 
HAPEX experiments and measurements are discussed by Andre et al. 
(1988).) 


Sensitivity tests of the stomatal resistance model 


Tests of the stomatal resistance model were done varying either soil water 
content, root-xylem resistance, or vapor pressure deficit. In response, leaf 
and epidermal water potential, stomatal resistance, and transpiration are cal- 
culated as functions of the imposed parameters. Solar flux was set =850 W 
m -2 , and vapor pressure deficit was set at 30 mbar. 


Soil water content 

Figure 4 shows the variation of bulk leaf and epidermal water potential, 
stomatal resistance and transpiration with soil water content. It is important 
to note that although the stomatal resistance eqn. (7) does not directly con- 
tain soil water content, stomatal resistance is still sensitively dependent on 
soil water content through iy t (see Appendix). Bulk leaf and epidermal water 
potential decrease with decreasing soil water content. Between 0. 1 2 and 0. 1 4, 
there are sharp decreases in leaf and epidermal water potentials as a result of 
a large decrease in y/ v and a large increase in Z g (a decrease in K t ). Stomatal 
resistance increases and transpiration decreases, however, where the epider- 
mal leaf water potential reaches the critical leaf water potential, which occurs 
at 9 V =0. 1 7. The value of d v at which critical leaf water potential is reached is 
dependent on the value of V. Higher values of V would cause the leaf water 
potential to fall below critical leaf water potential at even higher values of soil 
water content. Interestingly, transpiration appears to be almost independent 
of soil water content when y/ e is above y/ c . 
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Fig. 4. Soil water content (0 V ) vs. mesophyllic leaf potential (y/ x ), epidermal leaf water potential 
(^ c )» stomatal resistance (r # ), and transpiration (£*£*-). (Parameters and coefficients are the 
base values given in Table 4. ) 



z p 

Fig. 5. Root-xylem resistance (Z p ) vs. mesophyllic leaf potential (^), epidermal leaf water 
potential (^ e ), stomatal resistance (rj, and transpiration (L^f). (Parameters and coefficients 
are the base values given in Table 4. ) 
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Root-stem resistance 

The variation of bulk and epidermal leaf water potential, stomatal resis- 
tance and transpiration with root-stem resistance Z p is shown in Fig. 5. Bulk 
and epidermal leaf water potential decrease sharply with increasing Z p until* 
the critical leaf water potential is reached by y/ c at Z p equals 0.05. Above 
Z p =0.05, r s increases and transpiration decreases with increasing Z p . The large 
increase in stomatal resistance above Z p =0.05 prevents the water potential 
of the leaf from decreasing rapidly beyond this value of Z p . 


Vapor pressure deficit 

Figure 6 shows the variation of bulk and epidermal leaf water potential, 
stomatal resistance and transpiration with V. The most interesting result of 
this test is the corresponding changes in the calculated parameters when y/ e 
reaches the critical point at V= 32 mbar. As in the two previous Figures, sto- 
matal resistance increases sharply with further decrease in below y/ c ; tran- 
spiration, however, decreases with increasing V beyond y/ c because y/\ in- 
creases. The divergence of y/\ and y/ c beyond the critical leaf water potential is 
consistent with the feedforward observations of Farquhar (1978 ) and Sheriff 
(1984). In this model, the rate of divergence is dependent upon ft. It should 
be noted that the value of Fat which reaches y/ c is dependent upon (0 V ) 
and Z p . 



Fig. 6. Vapor pressure deficit ( V) vs. mesophyllic leaf potential epidermal leaf water po- 
tential (^ e )» stomatal resistance (r f ), and transpiration {LJE f ). (Parameters and coefficients 
are the base values given in Table 4. ) 
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TABLE 4 


Values of the coefficients, y/ c , Z p and M 


Variable 

Base 

Test 1 

Test 2 

Test 3 

Units 

^min 

210 

•* 

** 

** 

s m _1 

$ 

0.1 

0.0 

0.05 

0.1 

bar mbar" 1 

t >, 

-0.0001 

** 

** 

** 

s m -1 bar -1 

b 2 

o 

o 

** 

** 

** 

s m _l bar' 1 

C x 

0.5 

** 

** 

** 

sm' 1 (Wm- J )-' 

Cl 

10 

** 

*« 

** 

stir 1 (Wm- 2 )-' 

Vc 

-13 

-8 

-13 

-18 

bar 

Sc 

225 

*♦ 

** 

** 

Wm -2 

z p 

0.044 

0.03 

0.05 

0.07 

bar (Wm- 2 )' 1 

M 

0.65 

0.35 

0.50 

0.65 



Values of the coefficients used in the stomatal resistance (base) and boundary layer inclusive 
tests are shown in Table 4. Initialization is that of Lubbon, 16 June. 


Sensitivity tests with the boundary layer model 

The purpose of this section is to discuss what happens when the stomatal 
resistance model is implanted into the boundary layer model and the com- 
bined models are allowed to vary as a function of time for different initial 
conditions. More precisely, dependence of leaf water potential, stomatal re- 
sistance, transpiration, and radiative temperature on moisture availability, 
root— stem resistance, critical leaf water potential, and fi ( feedforward ) will 
be shown. The initial conditions are for HAPEX, 16 June 1986; the parame- 
ters and coefficients for these simulations are given in Table 4. 

Sensitivity to soil water content ... „ 

Changes in soil moisture availability and therefore in soil moisture affect 
the plant water flow eqn. ( 1 ) through Z, and y/ t , whose values change expo- 
nentially as 6 V decreases. It is apparent in Fig. 7a that there is little difference 
in the transpiration rate between M=0.50 and 0.65; at these moisture avail- 
abilities, y/ e is greater than y/ c . With a further decrease in y/ s , y/ e reaches 
leading to the formation of a transient transpiration plateau and then to a 
collapse of the transpiration as y/ t reaches y/ e . It should be noted that the col- 
lapse of the transpiration during the entire day is associated with an increase 
in the radiative temperature of 8K from that for M— 0.50 (Fig. 7b). 

Perhaps a useful means of representing water limitation is through WPSI. 
WPSI is much greater than zero for Af=0.35 (Fig. 7c); by comparison, when 
the supply of soil moisture to the leaves is sufficient to meet the demand, 
WPSI is less than zero. 
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T i ME (hr) 

Fig. 7. Sensitivity tests of (a) transpiration (Le£ f ), (b) radiative temperature, and (c) the 
Water Potential Stress Index (WPSI), varying the moisture availability ( M ). (Base coefficients 
are in Table 4.) 
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T I ME (hr) 


Fig. 7 contd. 


Sensitivity to Z p 

As shown in Fig. 5, larger values of the root-stem resistance, Z p (bar (W 
m -2 ) -1 ), result in a decrease in y/ x and y/ e . Figure 8a shows that even at high 
moisture availabilities, a relatively high value of Z p can cause a sharp increase 
in r s . This leads to a transpiration plateau which is about 100 W m -2 below 
that for a small value of Z p (Fig. 8b). In this case, the difference in transpir- 
ation corresponds to a canopy radiometric temperature difference of 2K. 
These results suggest the importance of Z p in plant models and its effect on 
boundary layer fluxes and radiometric temperatures. 

Sensitivity to y/ c 

Figure 9a, b and c shows how stomatal resistance, transpiration, and can- 
opy radiometric temperature vary with y/ c . The resistance is largest, and tran- 
spiration is smallest for — 8 bar. In comparison, for y/ c = — 1 8 bar, does 
not reach critical; thus r s resembles a U-shaped curve and there is no appear- 
ance of a transpiration plateau. This variation of r s with y/ e also leads to a 3K 
difference in simulated radiometric canopy temperature. 

The LPSI increases as increases and is >0 for y/ c ^ — 13 bar (Fig. 9d). 
Values of LPSI > 0 suggest a “leaf water deficit” and this deficit increases 
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7 9 1! 13 15 17 

TIME (hr) 


Fig. 8. Sensitivity tests of (a) stomatal resistance (r % ), and (b) transpiration (L,Ef), varying 
t he root-xy lem resistance ( Z p ) ( bar ( W m - 2 ) - 1 ) . ( Base coefficients are in Table 4. ) 
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with increasing y/ c . Note, however, that Fig. 9d implies that is larger for 
larger values of> c , suggesting that this deficit is relative to plant species and 
type. Perhaps, y/ c may be indicative of a mechanism by which the plant exer- 
cises direct control over the transpiration, thus preventing, in the short term, 
the leaf water potential from falling to a level that can harm the plant. 


Sensitivity to P , 

One other means of plant control of transpiration is through the feedtor- 
ward effect. Figure 1 Oa shows that transpiration decreases with £, and that 
there is a bi-modal transpiration curve for £=0.1 bar mbar -1 . This distribu- 
tion is associated with variations of y/ e , r s , and y/\ (Fig. 10b, c, d) that exem- 
plify results shown in Fig. 6. 

HAPEX measurements 

The HAPEX-MOBILHY (Hydrologic Atmospheric Pilot Experiment and 
Modelisation du Bilan Hydrique) program (Andre et al., 1986, 1988) was 
conducted over an agricultural region in the south of France during the sum- 
mer of 1 986. Its main purpose was to study the hydrological budget and evap- 
oration fluxes at the scale of a General Circulation Model. A network of rain 
gauges, flux-measuring devices and conventional meteorological observa- 
tions, supported by aircraft and additional ground surface measurements were 
made during an intensive phase which lasted from May until the early part o 

July 1986. ^ , 

One of the programs in HAPEX was directed jointly by Penn State and 

members of the Centre de Recherches en Physique de l’Environment (CRPE) 
and the Institut National de Recherches Agronomiques (INRA) to study the 
evolution of the stomatal resistance, leaf water potential and substrate water 
content during the intensive phase of HAPEX. Briefly, this experiment in- 
volved measuring the leaf water potential (using a pressure chamber), the 
leaf resistance to water vapor flow (using a porometer ) and the substrate water 
content with gypsum blocks, supplemented by gravimetric measurements. 
These measurements were made at two fields, one near the town of Lubbon 
and the other near Castelnau. The Lubbon field consisted of com, planted 
about the end of May; because of wet weather, the com at Castelnau was not 
planted until later in the field experiment. Periodically, LAI and plant heights 
were also determined. 

Surface energy flux measurements were made routinely at about a dozen 
sites by the Systeme Automatique de Mesure de l’Evapotranspiration Reelle 
(SAMER), including the Lubbon com field (Andre etal., 1988). Additional 
soil water measurements were made in various fields using neutron probe and 
gravimetric sampling; these are also discussed by Andr6 et al. ( 1 )• 

The importance of these HAPEX measurements with regard to this paper 
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ig. 9. Sensitivity tests of (a) stomatal resistance (r,), (b) transpiration (£*), (c) radiative 
mperature, and (d) the Leaf Potential Stress Index, varying the critical leaf potential (vO 
jar). (Base coefficients are in Table 4.) 
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Fig. 9 contd. 
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Fig. 10. Sensitivity tests of (a) transpiration {LJE { ) y (b) epidermal leaf water potential (^ c ), 
(c) stomatal resistance (r,), and (d) mesophyliic leaf water potential (w), varying^ (bar 
mbar -1 ). (Base coefficients are in Table 4.) 
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Fig. lOcontd. 
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TABLE 5 

Parameters for June 1 6 ( Lubbon ) 

Variable 

Lubbon 

Units 

^"min 

125 

s m~ 1 

/? 

0.02 

bar mbar" 1 

b, 

-0.0001 

s m“ 1 bar - * 

b 2 

500 

snr 1 bar -1 

c, 

0.15 

s m~ 1 (Wm- 2 )- 1 

C 2 

100.0 

snr 1 (W m~ 2 ) -1 


-13.0 

bar 

S c 

50 

Wm~ 2 

Zp 

0.035 

bar (Wm -2 )- 1 

Values of r min , and the coefficients used in running the stomatal resistance model for the HAPEX 
experiment. Values of critical leaf water potential and internal plant resistance (Z ) are also 
given. p 


is that they show the combined effects of evaporative demand, soil moisture, 
and root-stem resistance on stomatal resistance, leaf water potential, and 
transpiration. Furthermore, the data suggests the effect of root-stem resis- 
tance on stomatal resistance over a seasonal period. 

In simulating these results, values of the coefficients in eqns. (7), (8 ) and 
(9) were based loosely on measurements reported in the literature for corn 
(e.g. Choudhury, 1 983 ) of r % vs. solar flux and leaf water potential. The value 
of Z p was chosen on the basis of work performed by Abdul- Jabbar et al. 
(1984); see Table 5. 

Stomatal resistance 

Stomatal resistances (more properly, leaf resistances) were measured with 
a porometer during a period of about one month in HAPEX. Several leaves 
were sampled from different plants several times a day. Measurements were 
made on both adaxial and abaxial sides of the leaf and in three areas on each 
side, the tip, middle and stem region. Measurements were then linearly aver- 
aged over the several leaves for each of the leaf domains and then each leaf 
domain was combined into a single one-sided value of leaf resistance. The 
latter operation involved using different weights for the different leaf sectors, 
rather than taking a straight average of inverse resistances. (Mean values of 
r s are similar to straight inverse averages. ) 

All the stomatal (leaf) resistance data for the summer field program are 
shown in these two figures, one being for unstressed leaves (Fig. 11a) and the 
other for visibly stressed leaves (Fig. 1 lb), as determined subjectively in the 
field at the time of measurement. Note that the former exhibits a U-shaped 
distribution in the hourly averages, with a minimum at about 1 4:00 LST, while 
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Fig. 1 1. Normalized measured (relative) stomatal (leaf) resistances (small filled circles) for 
(a) unstressed and (b) stressed com leaves composited for the entire month of June 1986 (Field 
N6) at Lubbon, France. Circled crosses refer to hourly averages. Normalization refers to divi- 
sion by the minimum leaf resistance for the day. The average minimum stomatal resistance for 
June is indicated in the box. 
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Fig. 12. Simulated leaf resistance vs. time plotted against the measured leaf resistance (HAPEX) 
on 1 6 June 1 986, Lubbon. Measured resistances have been scaled by the shelter factor and LAI. 

the measurements for stressed leaves show a flat distribution, with a mini- 
mum at about 12:30 LST, and generally larger values of leaf resistance 
throughout the middle of the day than for non-stressed leaves. The difference 
between Figs. 1 la and b may be due to transient water stress. 

We also compared the measured values of stomata! resistances for 16 June 
in HAPEX. Figure 1 2 shows the measured resistances which have been scaled 
by LAI and the shelter factor P, vs. the simulated values, which are calculated 
for the leaves at the top of the canopy. The model appears capable of gener- 
ating a rich complexity of behavior, which is also found in actual measure- 
ments. For example, both simulated and measured resistances suggest a mid- 
day increase in leaf resistance, followed by a slight decrease. It is possible that 
the midday increase in the stomatal resistance is a manifestation of the leaf 
water potential reaching its critical value. 

The leaf water potential and transpiration plateau 

Much of the apparent increase in r , during the period of maximum solar 
flux occurs in conjunction with changes in leaf water potential. This is appar- 
ent in Fig. 1 3, which shows the variation of y/ x with time. Note that there is a 
period from about 12:00 to 18:00 LST when the differences between the leaf 
water potential and the predawn (base) value of the leaf water potential were 
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Fig. 13. Measured difference in leaf water potential from the pre-dawn (base) potential value 
(Jy/j) plotted vs. time. Dots represent individual HAPEX measurements composited for the 
month of June 1 986. Circled crosses represent hourly averages. 

almost constant with time in the hourly averages of the composite data. This 
behavior suggests that the leaf water potential was reaching its critical poten- 
tial during the afternoon, leading to a commensurate increase in r s . Since the 
base potential was typically between — 1 and - 2 bar, the leaf water potential 
distribution suggests a value of about — 1 3 bar for the critical leaf potential. 
Although we have defined in terms of a surface (epidermal) potential, the 
difference between the measured leaf (mesophyllic) and the epidermal leaf 
potential would not be large in these measurements, since the vapor pressure 
deficit effect in corn is thought to be small (B. Acock, personal communica- 
tion, 1989). Typically, corn seems to have a critical leaf water potential of 
about — 15 bar (B.J. Choudhury, personal communication, 1989). 

The appearance of a plateau in y/\ is also associated with a plateau in the 
transpiration (assumed to be proportional to the evapotranspiration). As seen 
in Fig. 1 4, both distributions of evapotranspiration measured by SAMER in 
adjacent fields were virtually identical before 08:00 h and after 1 6:00 h. In 
between, each field exhibited various degrees of the transpiration plateau. The 
dry field exhibits a wide plateau at about 300 W m -2 because both soil mois- 
ture is low and evaporative demand ( V) is high. A slight plateau may also be 
present in the wet field because of high evaporative demand; thus this plateau 
is of shorter duration and the magnitude of the flux is larger than in the unir- 
rigated field. 

The appearance of a leaf water potential and transpiration plateau is con- 
sistent with simulations shown in Figs. 8 and 9. If r s increases sharply at a 
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Fig. 14 Measured evapotranspiration (W m 2 ) during the HAPEX experiment, Castelnau, 
France. Dashed line is for a non-irrigated field of com while the solid line is for an adjacent 
irrigated field of com. Horizontal arrow indicates an evapotranspiration plateau near 300 
m -2 for the unirrigated field, suggesting that y t < during this period. 

critical value of leaf water potential, then a plateau in both the leaf water 
potential and transpiration must occur. If a transpiration plateau consistently 
occurs at a particular value of leaf water potential, then there must be a criti- 
cal value of leaf water potential owing to the constraints of the system (eqn. 

( 10) ). If one accepts that there is a critical value of leaf water potential, then 
there is no critical value of soil water potential (excluding a “nominal wilt- 
ing” value of soil water potential) or vapor pressure deficit at which r s in- 
creases sharply and transpiration levels off, because there is no monotonic 
relationship between leaf water potential and vapor pressure deficit or sou 
water potential. 

Variation of Z p . . 

We have plotted all the data for the minimum daily stomatal resistances 
(the values between 1 1:00 and 13:00 LST) and their hourly averages in a 
composite (Fig. 15). Regardless of how the choice of weighting functions for 
determining the mean r, was made, there was an unmistakable trend down- 
ward with time in the leaf resistance during the period, from about 265 to 1 60 

The smoothed decrease in the minimum daily leaf resistance shown in Fig. 
1 5 is intriguing because there was no correlation between the substrate water 
content and leaf resistance. Indeed, the water content in the root zone re- 
mained relatively high (about 0.18 cm 3 cm" 3 ) during June 1986. During this 
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Fig. 1 5. Average daily minimum stomatal (leaf) resistances (s m 1 ) for com leaves at approx- 
imately 1300 UT in a Field at Lubbon, France, during June 1986, as a function of date. Solid 
and open circles refer to alternate averages over leaf using the weighting factors shown above 
(values refer, from left to right, to the top, middle and bottom of the upfacing side of the leaf). 
The solid line is a subjectively determined fit through the data. 



Fig. 16. Simulated values of Z p (bar (W m~ 2 )- 1 ) vs. LAI and logarithm of the height of the 
com (m), based on Fitting the modeled stomatal resistance to the smoothed curve for measured 
resistances in Fig. 15. 


month, however, the plants grew from small shoots to over a meter in height 
and the LAI exceeded 3.0 in the corn by the end of the month. We suggest 
that the decline in r, with increasing LAI and crop height is related to the 
growth of prop roots, associated with an increase in the water transport effi- 
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ciency from soil to leaf. Such an increase in efficiency would correspond to a 

decrease in Z p . , , . . . r . 

This hypothesis was examined using the boundary layer model in the to - 

lowing wav. First, the smoothed trend of r s in Fig. 1 5 was plotted as a function 
of the LAI and crop height. Next, a value of Z p was chosen such that the 
simulation yielded the same daily minimum stomatal resistance measured on 
16 June at Lubbon. Using the same initial conditions, Z p was varied over a 
succession of simulations such that the minimum stomatal resistance gener- 
ated by the model followed the smoothed trend in minimum stomatal resis- 
tances depicted in Fig. 15. The result, shown in Fig. 16, is that Z p varied in- 
versely with both LAI and the logarithm of the crop height. This suggests the 
possibility that minimum daily stomatal resistance depends on the value ot 

z P . 


CONCLUSION 

The sharp increase in stomatal resistance associated with a decrease in ep- 
idermal leaf water potential beyond the critical leaf water potential can lead 
to the appearance of a transpiration plateau, or a bi-modal transpiration curve. 
Transpiration (and radiometric canopy temperature) is tied not just to the 
supply of moisture in the root zone and atmospheric demand, but to con- 
straints imposed by plant physiology which affects water flux through the 
plant. 


Suggestions for future research 


It is necessary to further examine the relationship between stomatal resis- 
tance, transpiration, and leaf water potential. There is also a need to clarity 
the vapor pressure deficit effect and to examine the role of the epidermal leal 
potential, and its relationship to the mesophyllic leaf potential Other factors, 
such as the effect of C0 2 and temperature on r s need to be further investi- 
gated, and a larger data base generated for all parameters. 

Equation (10) was derived under the assumption of steady state. Jones 
(1978, 1983), Waring and Running (1978), Running (1980b) and Tyree 
( 1988 ) suggest that the capacitance or storage capacity of the plant may be 
important. The effects of capacitance can be included in this model. 

Finally, the model is based upon a vertically averaged canopy with the im- 
plication of horizontal homogeneity. We question whether a vertical average 
properly captures the complex vertical variation of stomatal resistance within 
the canopy. More attention needs to be focused on the methods of making 
vertical and horizontal averages. 



STOMATAL RESISTANCE MODEL 


39 


ACKNOWLEDGEMENTS 

The authors would like to thank the sponsors of this research, the National 
Aeronautics and Space Administration under grant NAG 5-919, the United 
States Geological Survey under grant # 1 4-08-000 1-G 1490, and the United 
States Department of Agriculture/ Agriculture Research Service under grant 
#58-32U4-8-27. We are also grateful to Dennis Baldocchi, Basil Acock, James 
Bunce, Baskhar Choudhury, John Norman and Kenneth Shackel for invalu- 
able information and advice provided during the course of various discussions. 

APPENDIX: DERIVATION OF THE STEADY STATE EQUATION FOR 

Using eqn. (10) as the representation of steady state (Jarvis, 1976), we 
substitute eqn. ( 9 ) into eqn. ( 1 ) and eqns. ( 6 ) , ( 7 ) and ( 8 ) into eqn. ( 5 ) to 
derive the following: 

A^ e 2 + By/ e + C=0 (al) 

Since eqn. ( 8 ) is discontinuous, there are two parts to the solution to (a 1 ) . If 
<// g ^ y/ tc then the roots to this equation are 

A = -f(S)f( T)f(C0 2 )b ] (r cut + r af ) (a2) 

B=/(5)/(r)/(C0 2 )(r cut + r af ) 

[ - '■min + *1 Wt - 0V- H- Z, aV/ ( r cu , + r af )} ] ( a3 ) 


^cut ^af 

C=/(5)/(D/(C0 2 )(r cut +r af ) 

( 'mm ){V t - fiV-H- Z, c rV/ ( r cul + r af ) } ( a3 ) 

+ 'cut 'af( n't - P v ~ H ) - 'cut Z, aV 
if if/ i ^ i// tc then 

A=/(5)/(r)/(C0 2 )i 2 (r cul + r af ) (a4) 

B=/(5)/(7)/(C0 2 )(r cut + r af ) 

[-('min+^c+k^c) (a5 ) 

- - 0V- H- Zt aVj ( r cut + r af ) } ] 


^cut ^af 
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C=/(S)/m/(CO : )(r cut +r af ) 

(r mm +biy/ c +b 2 y/ c ) (a6) 

{i]/ t — f}V—H—Z K oV/ (r cut + r af )} 

"b ^ cui f af ( V^g ~~ — r cu\Z\OV 

where a is the density of vapor (air) times the latent heat of evaporation. 
(This solution is available on computer diskette from the authors. ) 
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ABSTRACT 


C and EM ' an<1 Sch f mugge > TJ > 199 °- Remote estimation of soil moisture availability 

and fractional vegetation cover for agricultural fields. Agric. For. Meteoroi. 52; 45-69. 

The purpose of this paper is to present a method for using remote measurements to estimate vene 

vegeTatmncanoDiJ^A 6116 ^ ^ ^ ““ and 5011 $urface water con,ents *>r partfal 

Zr!, ?, ' Pnma,T t0OlS are a boundar y layer model with vegetation and substrate 

Ponentii and two image products; the variation of surface radiometric temperature vs normalized 

^ n C : V , e8ela, T ,ndM (NDVI) ’ and ,he S,andard devia ‘ ion °f radiometric Surface te“ u « 
vs. radiometnc surface temperature. The method is based on determining: ( 1 ) asymptotic values of 
radiometric surface temperatures for sunlit bare soil and for dense sunlit vegetation and (2 ) a 

, ,p . "” ND y i “ ds " rf “'^ 

field e.ZLT.m™ > ” <l mid ' “ “*«, during ,he French HAPEX 


UNIKUDUCTION 


Models for estimating the surface turbulent energy fluxes and the soil mois- 
ure generally depend on a sensitivity of the surface radiometric temperature 
o soil water content. Over bare soil, variations in radiometric surface tem- 
f Tart- Ure te f d i t0 . corre,ated with variations in surface water content 

^ n f ta ’ 9?7;S f mU88e > I978 ; Jackson . 1982). Various models have 
been constructed to exploit the relationship between surface temperature and 
soil moisture. When used in conjunction with remote measurements of sur- 
face temperature such as determined from a satellite, these models yield es- 
timates of the surface moisture availability, the surface energy fluxes and the 
thermal inertia (Carlson and Boland, 1978; Carlson et al., 1981- Price 1 982- 
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Hatfield et al., 1984; RafTy and Becker, 1985; Taconet et al., 1986; Wetzel 
and Woodward, 1987; Lagouarde and Choisnel, 1990). Some of these models 
also take into account a layer of vegetation (Taconet et al., 1 986; Wetzel and 
Woodward, 1987; Lagouarde and Choisnel, 1990). 

The relationship between soil moisture and surface temperature is vastly 
more complex over vegetation than over bare soil. Over vegetation there is a 
considerable amount of temperature variability owing to the structure of the 
vegetation canopy and particularly to the amount of bare soil viewed by the 
radiometer and exposed to the direct solar beam. Leaves tend to be cooler 
than the exposed bare soil because the intercellular airspaces are nearly satu- 
rated with water vapor, which is drawn from a relatively deep soil layer, the 
root zone. Bare soil temperatures reflect the soil moisture only over the top 
one or two centimeters (Idso et al., 1975). 

Emitted surface radiance constitutes a blend of radiances emitted from 
either shaded or unshaded bare soil and vegetation. Thus, radiometric surface 
temperature variations over vegetated surfaces may be the result of variations 
in the amount of bare soil visible to the radiometer. The situation is compli- 
cated further by the problems introduced by canopy architecture, which in- 
cludes the variation in solar elevation angle and viewing angle of the 
radiometer. 

A significant step in the direction of modeling sparse or patchy vegetation 
cover was taken by Shuttleworth and Wallace (1985) who adapted the Pen- 
man-Monteith equation (Monteith, 1975) to account for energy partitioning 
between crop and soil. In order to circumvent the inconsistency of using a 
one-dimensional model to represent horizontally inhomogeneous surfaces, 
Shuttleworth and Wallace introduced the idea of two asymptotic temperature 
limits, one for bare soil and the other for vegetation. Their model was ex- 
pressed as separate bare soil and vegetation components with identical am- 
bient conditions above the crop. Latent heat fluxes from each of the compo- 
nents were combined using weighting factors, which were expressed in terms 
of combinations of soil and plant resistances determined from a knowledge 
of crop height and leaf area index (LAI). 

Shuttleworth and Wallace recognized that the intractable complexity of a 
vegetation canopy required that they make simplifications that allow them to 
avoid confronting the aspects of the three-dimensional canopy. While not dis- 
missing the effects of vegetation structure on the interception of solar radia- 
tion, they proposed that much of the fine-scale detail of the three-dimensional 
canopy would tend to cancel over time, leaving the first-order effects to be 
represented by a one-dimensional model. 

Stated differently, the fluxes over bare soil are relatively independent of the 
adjacent vegetation (and vice versa), but are linked via common substrate 
and atmospheric properties. In such a geometry, LAI has two values, one for 
the ensemble of vegetation clumps and one for the ensemble of vegetation 
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clumps and contiguous bare patches. The one-dimensional model operates 
separately in each regime; bare soil and vegetation. Total fluxes for the com- 
bined bare soil and vegetation are related through an additional parameter, 

the fractional vegetation cover. ... 

We propose a similar approach to that of Shuttleworth and Wallace, in which 
a boundary layer model is used in conjunction with remote measurements ot 
surface temperature and normalized difference vegetation index (NDV1) to 
calculate fractional vegetation cover, LAI, substrate water content in two lay- 
ers and surface energy fluxes. The method is illustrated using remote mea- 
surements made during the French HAP EX/MOBILHY field experiment. 


THE BOUNDARY LAYER MODEL 

Like Shuttleworth and Wallace, we will ignore the intractable aspects of the 
three-dimensional canopy, and imagine a simple structure in which clumps 
of vegetation (which may consist of a single plant) are interspersed with bare 
soil patches (Fig. 1 ). Within the vegetation clumps, no direct solar radiation 
reaches the ground and no bare soil is visible to the radiometer. Uniformly 
distributed sparse vegetation and small dense clumps of vegetation are sepa- 
rated, perhaps randomly, by bare soil. The fractional part of the surface cov- 
ered by the vegetation clumps is/ v and that covered by bare soil is ( 1-A). 
No distinction is made between sunlit and shaded vegetation or between sun- 
lit and shaded soil visible to the radiometer, although shaded bare soil may 
be much cooler than unshaded bare soil. Vegetation fraction also depends on 
the viewing angle of the radiometer, it will effectively increase with decreas- 



Fig. 1. Schematic view of a partial vegetation cover with leaf area index (LAI) equal toLOm 
contiguous vegetation patches (hatching). Area represents that of a single pixel or combination 

of pixels. 
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ing elevation angle of the radiometer. This variation of f v with viewing angle 

is not considered. . • 

Since the details of the boundary layer model are not as important as the 
technique of applying such a model to the estimation of fractional vegetation 
cover and soil moisture, we will present only the basic elements of the archi- 
tecture and equations; the reader is referred to Carlson et al. ( 1 98 1 ) and Ta- 
conet et al. (1986) for a discussion of the bare soil and vegetation compo- 
nents, respectively. Results presented in this paper should be reproducible 
with any comparable model, e.g. that of Wetzel and Woodward (1987). 
Henceforth, the boundary layer model will be referred to as the CM. An ov- 
erview of the plant and substrate structure in the CM is presented in Fig. 2. 

The CM contains a single layer of vegetation, and interleaf air layer, two 
soil layers, a surface and a mixing layer. Fractional vegetation cover (/ v ) is 
specified for a mixture of bare soil and vegetation. In the vegetation compo- 
nent, vegetation density is expressed in terms of LAI. Bare soil and vegetation 
regimes are calculated separately, but allowed to interact through exchanges 
of momentum, heat and water vapor with a common surface air layer above 
the canopy and a common substrate below. 

Downward solar flux ( S ) and downward long-wave radiation (R,l) are 
identical over both bare soil and vegetation regimes. Surface fluxes of sensible 
and latent heat and upward fluxes of long-wave radiation for bare soil and 
vegetation (R lb T and R |V T ) are determined separately for each regime. Radio- 
metric temperature of the canopy ( T c ) is computed from a weighted average 
of bare soil and vegetation components of upward long- wave radiation fluxes. 

The flow of water vapor and sensible heat is expressed as follows. In the 
vegetation, latent heat flux (LE { ) passes from the inside of the leaf at leaf 
temperature ( T { ) and specific humidity (<?i(7i)) to the interleaf air spaces, 
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Fig. 2. Overview of plant canopy model. Fractional vegetation cover is/ v and the weighting for 
fluxes ( F) is given by the formula at the top. The remaining symbols are defined in text. 
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which have a specific humidity (q a f ), through the stomatal and the leaf 
boundary layer resistances (r s ; r af ). Stomatal resistance is calculated from a 
knowledge of substrate water content, incident solar flux and LAI. Sensible 
heat flux (H { ) passes from the leaf (at temperature ( 7", ) ) to the interleaf air 
spaces, which are at temperature ( T af ), through the leaf boundary layer resis- 
tance (r af ). 

At the soil surface below the plants, latent heat flux ( LE t ) passes from the 
ground surface beneath the canopy, which is at temperature ( T gv ). to the in- 
terleaf air space through surface resistance (r^) and interleaf air space resis- 
tance (r ag ). Sensible heat flux ( H % ) passes from the ground across resistors 
(r ag ) between the ground temperature and the interleaf air space temperature 
(T af ). These parallel fluxes (L£ f ; LE t ; H f ; H t ) sum to the total vegetation 
fluxes of latent and sensible heat (LE V \ H v ) at the top of the vegetation can- 
opy and pass through the canopy boundary layer resistance (r bf ) and the sur- 
face layer resistance (r a ). 

Similarly, the calculation of fluxes over bare soil patches is simpler, but 
expressed in the same notation as over vegetation. Latent heat flux passes 
from the substrate through a soil resistance (r g ), a boundary layer resistance 
(r b ) and a surface layer resistance (r a ). Sensible heat flux passes from the 
ground surface at temperature T gb directly to the atmospheric surface layer 
through the resistors (r b ) and (r a ). The reader is referred to Taconet et al. 

( 1 986 ) for further details concerning the model. 

For a mixture of solid vegetation and bare soil patches, fluxes of sensible 
and latent heat and the upward flux of long-wave radiation above the plant 
canopy and the substrate heat flux ( G ) are taken as weighted averages of the 
bare soil and vegetation components, according to the vegetation fraction (J v ) 
using the weighting equation shown in Fig. 2. Temperature and specific hu- 
midity at the top of the surface layer (T a , q a ) and the temperature and water 
content in the soil are identical for both bare soil and vegetation fractions. Air 
resistances for both vegetation and bare soil fractions depend on the wind 
speed and temperature above the canopy and therefore, indirectly, on condi- 
tions in both the vegetation and bare soil components. Radiometric temper- 
ature of the canopy (T c ) is determined from the weighted upwelling long- 
wave radiation flux. 

An important parameter in the model is the moisture availability, which is 
defined for the soil surface layer with volumetric water content d vo as 

M 0 = 0 VO / 0 sa . = ( r a + r b ) / ( r g + r a + r b ) 

= (r a + r br )/(r bg + r }l + r b f) (1) 

and for the root zone water content as 

M=6 v /d sa, 


( 2 ) 
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where 0^, is the volumetric water content of the soil at field capacity (0 V ). 
Moisture availability, the primary moisture parameter, is also defined for the 
surface air layer as the ratio of evaporation to potential evaporation. It should 
be noted that M a , is identical for both bare soil and vegetation dominants. 
Since r b and r bf are calculated, eqn. ( 1 ) imposes small differences in the val- 
ues of r g and r bg , which are not directly employed. Further, there is some con- 
tention in the idea that M a can be directly equated to soil water content, al- 
though it is generally agreed that the former depends closely on the latter. 

MEASUREMENTS 

The HAPEX field experiment 

The HAPEX-MOBILITY (Hydrologic Atmospheric Pilot Experiment and 
Modelisation du Bilan Hydrique) program was conducted over a 100 km by 
1 00 km square, located within an agricultural region over southwestern France 
during the summer of 1986 (Andre’ et al., 1988). Its main purpose was to 
study the hydrological budget and evaporation fluxes at the scale of a General 
Circulation Model. At network of rain gauges, flux-measuring devices and 
conventional meteorological observations, supported by aircraft and addi- 
tional ground surface measurements were made during an intensive phase 
which lasted from May until mid- July 1 986. 

One of the programs in HAPEX was directed jointly by Penn State (PSU), 
the Centre de Recherches en Physique de l’Environment (CRPE) and the 
Institut National de Recherches Agronomique (INRA) to study the evolu- 
tion of various plant parameters (stomatal resistance, leaf water potential, 
LAI, crop height) and their relationship to substrate water content and radio- 
metric surface temperature during the intensive phase of HAPEX. PSU ground 
measurements were made in two fields of com, one near the town of Lubbon 
(in Field N6) in the north of the test area and the other near the town of 
Castelnau in the south. The former was situated in a field of com with sandy 
soil, within a large clearing surrounded by a forest (Les Landes); the latter 
was located in a com field on a ridge. Soil water content was measured once 
or twice daily using gypsum blocks. These devices were calibrated at both 
sites by gravimetric methods, both in situ and in separate pots of soil. Addi- 
tional soil water measurements were made periodically by PSU, INRA and 
NASA groups at this and other sites using gravimetric sampling or neutron 
probe. Unfortunately, because of late planting due to wet weather, the com at 
Castelnau did not emerge until mid-June. Consequently, the focus of this pa- 
per is on the Lubbon measurements. 

The aircraft program consisted of a series of 16 flights with a NASA C-l 30 
aircraft over the target area during the intensive phase of HAPEX. Two multi- 
spectral scanners covering the wavelength region from the visible to the ther- 
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mal infrared and 21 -cm microwave radiometer were mounted aboard this 
aircraft. Fifteen flights were made at 300 m and 1500 m; additional high al- 
titude (6000 m ) flights were made for the purpose of comparison with satel- 
lite measurements. One of the scanners was the Thermal Infrared Multi-Spec- 
tral Scanner (TIMS) and the other one the NS001 scanner which was 
functioning at solar wavelengths only (Schmugge and Janssen, 1988 ). Overf- 
lights at Lubbon were made at approximately 1 1 : 30 h local time (LST). The 
focus of this paper is on 4 days during the HAPEX experiment, 6, 16 and 27 
June and 2 July. These were the only days in which flights were made during 
clear or nearly clear sky conditions. 

Radiometric temperatures 

Temperature measurements were obtained in an across-track scanning mode 
in the form of digital (DN) counts interpolated between cold and hot refer- 
ence sources. Absolute accuracy of the blackbody references are estimated to 
be about 0.5 K, resolution accuracy about 0.35 K and calibration accuracy 
about 0.8 K. Corrections for water vapor and carbon dioxide attenuation as 
functions of viewing angles were made using a model of Price ( 1982 ), which 
calculates the temperature correction for water vapor and carbon dioxide at- 
tenuation for different viewing angles given a temperature and moisture 
sounding. Since all the flight data discussed in this paper refer to the 1 500 m 
altitude passes, the total water vapor correction was only a few degrees K. 
Overall, the accuracy of the radiometer including atmospheric correction is 
felt to be about 1 .0- 1 .5 K. Pixel size for the 1 500 m flights was approximately 
4 m in diameter at nadir; no corrections for pixel size as a function of viewing 
angle were made. 

NDVI 

The NS001 Thematic Mapper simulator (TMS) was used to measure short- 
wave spectral radiances from the C- 1 30. Although having a similar resolution 
to that of the TIMS and an across-track scanning pattern, the instantaneous 
field of view and surface co-location was slightly different from the TIMS 
radiometer owing to different scan angles and swath widths. NDVI measure- 
ments were made using the formula: 

NDVI=(^4-^3)/(^4 + J?3) (3) 

where is the radiance from Channel 3 (wavelength band 0.633-0.697 pm 
and R a is that for Channel 4 (wavelength 0.767-0.910 pm). No corrections 
were made for solar elevation and radiometer viewing angles. Viewing angles, 
however, were all less than 30° . According to Paltridge and Mitchell ( 1 989 ) , 
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corrections to a nadir NDVI for NOAA satellites are less than 10% for view- 
ing angles less than 30°. 

% 

Soil moisture 

Eight gypsum block locations were installed at each of the two PSU field 
sites, Lubbon and Castelnau. Average and standard deviations of volumetric 
soil water contents ( 0 V ) were determined for each site at approximately 5,10, 
20 and 40 cm depth. In general, the standard deviation of 0 V was about ± 0.04. 
Although this is not an insignificant variation, it is very similar to the mean 
standard deviation obtained by Bell et al. (1980) for that of natural soils. 

ANALYSES 

General approach 

Fractional vegetation cover and LAI are determined by matching simu- 
lated radiometric surface temperatures, vegetation fraction and LAI with 
measured values of NDVI and radiometric surface temperatures. Asymptotic 
limits for sunlit leaves and sunlit bare soil are determined by inspection of 
the distribution of NDVI vs. radiometric surface temperature. Asymptotic 
leaf and soil temperatures are used to solve for, respectively, the root zone 
and soil surface water contents and fractional vegetation cover using the CM. 
Spurious effects of mutual leaf shading and small-scale variations in soil 
moisture are largely ignored, although they can be assessed qualitatively from 
inspection of image products. We make no distinction between different types 
of crops (com, soybeans) in either the modeling or the interpretation of the 
images. 

Method 

The significance of ND VI 

The NDVI constitutes a very useful tool in the study of vegetation cover. 
Normalization of the radiance differences between spectral intervals on either 
side of the near-infrared reflectance discontinuity tends to reduce effects of 
changing sun or viewing angle. Importantly, the NDVI is known to be corre- 
lated with the green leaf area, such as the one-sided LAI (Tucker, 1979; Hol- 
ben et al., 1 980; Curran, 1 983; Asrar et al., 1 984; Best and Harlan, 1 985; Gallo 
et al., 1985; Sellers, 1985; Hansen and Soegaard, 1987; Peterson et al., 1987; 
Nemani and Running, 1989a). Tucker (1979) and later Best and Harlan 
(1985) show that the optimum sensitivity to biomass and LAI can be found 
in the ratio of near infrared (0.70-0.80 pm) to red (0.63-0.69 pm) radi- 
ances, bands which are similar to those referred to above. When these ratios 
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are compared with actual field measurements the error in LAI is typically 

about ± 1.0 for a given site. inrrpasinc LAI and then enters an 

NDVI increases almost very slotly with increasing LAI. 

asymptotic regime in which N asymptotic region pertains to a sur- 

Curran ( 1983) points out Although there is some variation, 

face almost completely co^ere^by ^ density or LAI usually occurs 
an upper asymptote of ND 1 • S limit however is rather van- 

near 0.5-0.8 for dense vegetation. (Pal, ridge 

able and depends on vegetation type^age, between -0.1 and 0.2. 

and Barber, 1 988 ) . For bare soil ND VI t ends tc > ™ at values of 3- 
Curran also shows that the asymptotic region grasses. Asymp- 

4 for short crops such as wheat corn, sor ^J* Holben et al. (1980) for 

above which NDVI rapidly approaches an asymptotic 1 1 

however, the asymptotic region is in ex- 

bers of clumped needles. 0R4 , di<5CUSS contributions to the 

asymptotic part of the LAI ys. rsu increase occurs in the fraction of 

is close to 100%. Above an LAI of 3C L Mfle ..ncreasec . ^ ^ ^ ^ 

bare soil visible to the radiometer. A ®. (f = \ 0) begins at 

vegetation fraction (4) is that a*. d tmn ove ^ ^ of 

LAI = 3.0. LAI can increase above A nsofar as (he sensitiv ity of 

fu" A^ £ Se mS't be more appropriate for other types of vege- 
tation, such as trees. 

The vertical variation of soil moisture “iLubbrn nle Qf ^ 

Figure 3 is * ‘ at' considerable scatter, it is clear 

water content for 16 June at Lubbon. P re i a tivelv moist (0 V =O.2), 
that the deeper substrate (e ow^ ou should be noted that all of 

while the surface is very dry (0 v =O.UZ-u.uo;. 
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Fig. 3. Vertical profile of volumetric soil water content (0J vs. depth in cm measured by differ- 
ent devices on 16 June 1986 at Lubbon. 



Fig. 4. Variations in measured volumetric water content of soil (0 V ) at various substrate levels 
during June 1986 at Lubbon. The 5, 10, 25 and 40 cm curves represent measurements made 
with gypsum blocks in the N6 corn field. The 0-5 cm curves are from gravimetric samples in 
the oat (N2), com (N6)and bare (Nl) fields at Lubbon. 
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the gravimetric measurements show a marked soil surface dryness after the 
early part of June until after 23 June (Fig. 4). By contrast, the deeper sub- 
strate remained fairly moist, but with some decrease in the 5 and 10 cmr val- 
ues of d v between 10 and 19 of June. Some irrigation was carried out on and 
after 23 June owing to lack of precipitation during the first three weeks of the 
month. Despite the irrigation and some minor rainfall episodes during the 
last few days of June, the large vertical soil water content gradient was present 
as late as 27 June. By early July, surface soil water contents had increased 
toward those at deeper levels. 

Thermal imagery 

Figure 5 is a thermal image made from digital scan data from the C-130 
aircraft over the Lubbon region. Here, temperature is represented by a gray 
scale shading, dark representing cooler temperatures and light warmer tem- 
peratures. Field N6, consisting of corn about 0.8 m high with an LAI of 1.8, 
was the location of the Penn State measurement operation. The average tem- 
perature of this field was about 47 °C. In contrast. Field N2 contained oats 
with a temperature of about 28 °C, while Field N4 with newly planted corn 
and Field N 1 , which was bare, exhibit radiometric surface temperatures well 
in excess of 50 °C. This image clearly shows the large differences in radiomet- 
ric surface temperatures between the exposed sunlit bare soil areas and the 



Fig. 5. TIMS thermal infrared image made from NASA C- 1 30 flying at 1500 m on 16 June 1986 
near Lubbon. Dark is cool and light is warm. Field numbers are indicated and field averages of 
temperature ( °C) and NDVI are listed. 
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vegetation, which is uniformly dark. It indirectly reflects the large vertical 
variation in soil dryness between the surface and substrate shown in Figs. 3- 

4. 

Large temperature differences between vegetation and bare soil occurred 
on all 4 days in which temperature images were made. While less bare soil is 
visible on 27 June (figure not shown) than on 16 June (Fig. 5), the distri- 
bution of warm bare soil and cool vegetation temperatures are similar. By 
contrast, however, images for 6 June and 2 July (not shown) respectively 
exhibit small and large fractions of vegetation. 

Combining ND VI and thermal measurements 

Spatial variations in surface radiometric temperature, such as shown in Fig. 

5, are related to variations in the vertical variation of soil water content mod- 
ulated by fractional vegetation cover. The relationship between radiometric 
surface temperature and vegetation amount is illustrated in Figs. 6-9. Sym- 
bols plotted on the figures pertain to individual resampled pixels. The best 
example is that of 1 6 June, which is shown in Fig. 6. Data fall into two general 
clusters. One cluster contains positive values of NDVI in an envelope that 
slopes from the upper left-hand to lower right-hand sides of the NDVI/r c 
distribution (Fig. 6a). The other cluster comprises points scattering horizon- 
tally and having negative values of NDVI. The sloping region extends from 
cool temperatures and high NDVI over the oat field (Field N2) to very high 
temperatures and low values of NDVI over newly planted corn fields (Field 
NO). 

The figures exhibit a remarkable range of surface temperatures, with very 
warm values over bare or nearly bare soil. In the sloping distribution, pixels 
with high NDVI (in excess of 0.6) and low temperatures represent the 
asymptotic region in which the temperature approaches that for a sunlit leaf; 
(shaded leaves would be somewhat cooler but with approximately the same 
NDVI). At this end of the distribution, the upwelling long-wave radiation 
from the soil is almost totally attenuated. Large NDVI values for the oats 
probably correspond to an LAI well in excess of 3.0 and a vegetation fraction 
of 1 .0. At the other end of the distribution, the very warm temperatures cor- 
respond to virtually bare soil with LAI almost equal to 0 and a zero vegetation 
fraction. 

Fig. 6a. Distribution of NDVI vs. radiometric surface temperature ( C C, corrected for atmos- 
pheric attenuation ) for individual pixels over the scene in Fig. 5. The type of vegetation planted 
in each Field is indicated by a symbol defined below the figure. 

Fig. 6b. Same as Fig. 6a but in schematic form to show how the axis of variation (dashed line) 
was constructed from the data in Fig. 6a. Circled numbers are the values of LAI determined 
with the aid of the boundary layer model by matching observed and simulated temperatures 
given the surface and root-zone moisture contents. Decimal number represents the fractional 
vegetation cover, circled numbers the derived LAI, and the underlined number the value of 
surface roughness (cm ) used in the model for that simulation. 
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27 June 1936 Lubbon 


NDVI vs Tims Ts 



20 30 +0 50 60 


Tims Ch5 Ts (C, corrected) 

□ N1 5oH + N2 Oats O N3 Corn A N6 Corn 

Fig. 7. Same as Fig. 6a but for 27 June. 

Two other features in Figs. 6-9 are worth noting: one is a horizontal scatter 
of surface temperature, a cluster of pixels with a constant NDVI (e.g. the 
horizontal width of the envelope of data points in Fig. 6a). This cluster is 
probably associated with spatial variations in soil moisture and to variable 
effects of the canopy architecture, e.g. mutual shading of leaves. The second 
feature is the apparent decrease in the maximum NDVI with time between 6 
June and 2 July. This trend may be due to changes in vegetation color with 
age. 

Goward et al. ( 1985) and Nemani and Running ( 1989b) have proposed 
that NDVI/Tc diagrams can be used to infer vegetation amount or canopy 
resistance. The axis of the sloping distribution of points in Figs. 6-9 defines 
the variation of vegetation fraction with surface radiometric temperature. 
Providing that there is not a great deal of noise and the distribution of points 
exhibits a definable slope over a wide range of NDVI, one can define an “axis 
of variation” through the center of the sloping distribution of points, which, 
in this case, connects the mean values of T c and NDVI for each field. This 
axis describes a spectrum of temperatures and NDVI values with visible bare 
soil fractions varying from almost zero to virtually 100% cover. For the sake 
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NDVI vs Ts 


July 2 lubbon Hapex 



Tims Ch5 Corrected Surface Temperature 
□ Field 2 Oots + neld 6 Com 

Fig. 8. Same as Fig. 6a but for 2 July. 


of illustration, the axis of variation for Fig. 6a is shown by itself in Fig. 6b, 
but superimposed on a schematic background of pixel represented by hatching. 

The NDVI/Fc distribution was less coherent on 27 June (Fig. 7) than on 
16 June, although the variation in surface temperature was as large. A gap 
exists in the axis of variation at high surface temperatures and low values of 
NDVI owing to the relatively few bare soil patches between the plants. The 
separate cluster of bare soil temperatures, representing bare fields, is still 
present. 

Field N6 was relatively cool and verdant on 2 July. Considerable scatter in 
the ND VI /T distribution (Fig. 8) may have been caused by the effects of 
irrigation, which was administered unevenly, and to ripening or senescence 
of the oats. Consequently, no axis of variation could be drawn with much 
confidence. 

Most of the fields on 6 June were bare. The NDVI/T C distribution (Fig. 9 ) 
for that date shows a horizontal variation of radiometric surface temperature 
along the zero NDVI axis, with one isolated cluster of pixels at NDVI =0.6 
and 10°C. No axis of variation could be drawn with any confidence. 

An oddity of Fig. 6 is the presence of a separate horizontally distributed 
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Fig. 9. Same as Fig. 6a but for 6 June. 

cluster of pixels near or below the zero NDVI axis. This cluster represents 
pixels from unplanted or newly planted fields without vegetation, as distinct 
from bare soil patches between plants in fields that contain maturing vegeta- 
tion. The reason for this difference is that freshly plowed bare fields would 
exhibit clods of relatively wetter earth interspersed with fragments of very dry 
surface soil, whereas bare soil between the plants would be dryer and littered 
with dead plant debris. 

Interpretation of the imagery using a boundary layer model 
Based on the assumptions described above, fractional vegetation as a func- 
tion of NDVI or surface temperature was derived with the aid of the bound- 
ary layer model according to the following procedure. 

First, asymptotic temperature extrema were determined for sunlit bare soil 
and sunlit leaf temperatures with the aid of the NDVI/7’ C distribution (Figs. 
6-9) and the “arch diagram”. The arch diagram was originally developed by 
Coakley and Bretherton ( 1982) to analyze cloud cover in partially filled fields 
of view and later used by Albrecht et al. (1988) to separate sea surface tem- 
peratures from those at the top of marine stratus clouds. Figures 10-12 show 
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Fig. 10. The arch diagram: standard deviation of radiometric surface temperature vs. uncor- 
rected radiometric surface temperature ( °C ) for 16 June at Lubbon (the image in Fig. 5 ). Each 
data point represents a sub-area consisting of about 10X 10 original pixels. 


the arch diagrams for 16 and 27 June and 2 July. The arch diagram consists 
of the standard deviation of temperature (in 10x10 pixel subsets) of the 
image vs. radiometric surface temperature averaged over that subset. No cor- 
rection was made for viewing angle or for atmospheric attenuation. The two 
“feet” of the arch, characterized by dense clusters of relatively low standard 
deviations, suggest areas of uniform vegetation and bare soil temperatures, 
which correspond to the temperature extrema along the axis of variation in 
the NDVI/T C distribution (e.g. Fig. 6b). 

The second step was to use the boundary layer model to calculate substrate 
water content using the asymptotic temperatures. Thus the water contents 
(0 VO ; d v ), pertain, respectively, to sunlit bare soil (f v = 0; LAI = 0; T Q =55°C) 
and sunlit vegetation (/ v = 1.0; T c =28°C). Solutions are obtained by forcing 
agreement between simulated and measured surface temperatures at flight 
time (11:30 LST) (Table 1 ) for bare soil and solid vegetation. 

The third step was to use these derived substrate water contents to calculate 
/ v as a function of surface radiometric temperature. In these simulations, LAI 
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Fig. 1 1 . Same as Fig. 1 0 but for 27 June. 


is fixed at 3.0 for the vegetation fraction and 0 for the bare soil fraction (Fig. 
1 ). We force the simulated and measured radiometric surface temperatures 
to agree along the axis of variation for different values of / in the model. 
Results are not sensitive to the choice of the solid vegetation value for LAI 
since the model is used only to convert a given temperature range to a fraction 
(/ v ) between 0 and 1 .0. Given / v vs. T cy one obtains the distribution of NDVI 
vs. / from the NDVI/T C diagram. Figure 6b shows the values of LAI, / and 
surface roughness used to simulate temperatures along the axis of variation 
on 16 June. For example, given LAI = 1,8 (Field N6) the radiometric surface 
temperature is about 46 °C, the NDVI is about 0.4 and/ = 0.6. 

A final step was to perform a simulation to generate a range of surface ra- 
diometric temperatures as LAI varies above 3.0 with /= 1.0. This operation 
simply allows us to determine LAI as a function of radiometric temperature 
for the solid cover, but has no bearing in determining/,. Similar calculations 
were made for 27 June. 

At this point, it is possible to specify vegetation fraction (/ between 0 and 
1.0) and the average value of LAI for each resampled pixel, given its mea- 
sured values of radiometric surface temperature and NDVI. From this infor- 
mation one calculates the fluxes of sensible and latent heat for each pixel. 
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Fig. 1 2. Same as Fig. 10 but for 2 July. 


TABLE 1 

Average values of measured and deri ved quantities for Field N HAPEX (Lubbon), 1986 

“><">• * [4r f”> k> NDV,,n ” "■ 


6 June b 
16 June 
27 June 
2 July 


0.6 46 

1.0 33 

1.0 32 


; " measured; -No aircraft measurements made over Field N6; 'Values in parentheses are uncenam. 


DISCUSSION 

Results of simulations 

Not surprisingly, the high temperature asymptote corresponds to a low soil 
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surface water content (Af 0 = 0.02; 0 V = 0.006). We were unable to simulate the 
highest temperatures without specifying a very low thermal inertia, appropri- 
ate to the sandy soil of Lubbon. Simulations for wet soil show little variation 
in temperature as a function of M between M= 0.6 and M— 1.0 (6 V between 
0.2 1 and 0.34). Thus, it was more difficult to simulate precisely the root zone 
water content at the cool end of the axis of variation, since surface tempera- 
ture is insensitive to substrate water content in wet conditions. Accordingly, 
we chose to let A/=0.65 and 0 V = 0.22 for both 16 and 27 June, values corre- 
sponding to those measured in the root zone (Figs. 3-4). 

Simulations were made for the bare soil cluster (NDVI less than zero). 
These showed that the surface moisture availability ( M 0 ) varied from about 
0 to 0.6 between the warm and cool temperatures. 

Asymptotic temperatures on 27 June were approximately the same as on 
1 6 June, although it was more difficult to derive a bare soil asymptote on 27 
June because of greater scatter and a gap in the measured NDVI/T C distri- 
bution. The arch diagram (Fig. 11) for 27 June suggests a small amount of 
dry, bare soil with temperatures approaching 60°C, but the NDVI/T C distri- 
bution does not lend confidence in determining the warm end of the axis of 
variation. (Recent evidence published by the French indicates that the mea- 
sured temperatures from TIMS are too warm, perhaps by as much as several 
°C for the high temperatures. This error is due to deterioration in the block 
body target in the radiometer. See: M. Stoll, Temperature de Surface restituee 
a partir du radiometre aeroporte T.I.M.S. Rapport de stage effective au 
CNRM, Toulouse, 29 June 1989 promotion I.M. 87/89, 56 pp.) Simulations 
for 27 June yield similar results to those of 1 6 June, which is not surprising in 
view of the fact that the NDVI/T C distributions appear to be similar. How- 
ever, we cannot explain why the axis of variation was at lower values of NDVI 
on 27 than on 16 June (Fig. 7). 

Simulations for 6 June show that the surface moisture availability ( M 0 ) 
varied from about 1.0 on the cool side to about 0.1 (0 V =O.O3) on the warm 
side of this cluster of bare soil temperatures. Surface temperatures were not 
as warm as during late June, although a few pixels were extremely hot. The 
isolated cluster of cool pixels in Fig. 9, which represent the oats canopy, could 
not be simulated without reducing the solar flux intensity below that for clear 
skies. This suggests that these pixels may have been partly shaded by leaves 
or by cloud. 

Although no axis of variation was obtainable for 2 July because of the enor- 
mous amount of scatter, the NDVI/r c distribution does suggest that the 
warmest soil temperatures were about 45 °C, the equivalent in the model of a 
surface moisture availability ( M a ) of about 0.1. These values reflect the 
somewhat higher surface soil moisture contents found on this day (as com- 
pared with 16 June); the arch diagram indicates that there were still some 
very high surface temperatures corresponding to almost zero surface moisture 
availability in some bare patches. 
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Table 1 summarizes some results for Field N6. Presented are values of LAI 
(measured), LAI (derived),/, (derived), NDVI (measured) and the pure 
sunlit leaf and sunlit bare soil temperatures (LAI = 3.0 and 0) and for the 
appropriate field average of LAI which allows measured and simulated sur- 
face temperature to agree. Also shown are simulated root zone water contents. 
The last two days constitute 100% canopy cover for Field N6 (LAI exceeding 
3.0), whereas 16 June is an example of partial vegetation cover in a single 
field. 


Comparison with measurements 

Some independent confirmation of the results can be made by comparing 
the measured values of LAI and NDVI with the simulated values of NDVI 
and LAI as a function of simulated vegetation fraction (/ v ). LAI was mea- 
sured at the PSU site in Field N6 on four occasions during the same period. 
Figure 1 3 shows the distribution of measured NDVI vs. measured LAI and 
the derived NDVI (from the axes of variation) as simulated by the model. 
The comparison is obviously very good for 16 June but agreement with mea- 
surements is not as close for 27 June because of greater scatter in the data 
(Fig. 7). Had we chosen a somewhat smaller LAI threshold for solid vegeta- 
tion than 3.0, agreement between measured and modelled LAI would have 
been closer to the measured LAI on 27 June. 



Fig. 13. NDVI vs. percent cover and LAI for model simulations for 16 June and 27 June at 
Lubbon. Squares represent simultaneous aircraft measurements of NDVI and ground measure- 
ments of NDVI made on 4 days over the Lubbon field (N6) during June; crosses represent 
measured NDVI and LAI for almost bare soil. Heavy solid line represents simulated LAI and 
fraction vegetation cover vs. NDVI, consistent with axis of variation for measurements. The 
thin solid line obeys a formula proposed by Gallo et al. ( 1 985 ) for NDVI vs. LAI over com. 
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Problems in applying the method 

% 

The difficulty in determining the axis of variation, especially on 6 June and 
2 July underscores the central problem in this method. To some extent, poorer 
agreement with measurements on 27 June, as compared with 16 June, is due 
to the wide horizontal scatter of points and the lack of values on the warm 
side of the axis of variation. Clearly, the distribution of points on 16 June is 
ideal, because the variance of temperature at constant NDVI is quite small. 

The derived fractional vegetation cover depends on the architecture of veg- 
etation, the geometry of bare soil distribution, the type of soil and vegetation 
and the solar and radiometer viewing angles; will tend to increase with in- 
creasing solar or viewing angles. The present results, however, suggest that 
various crops can be combined to produce a meaningful axis of variation. 
However, extending the method to a mixture of trees and crops may prove 
untenable. 

A difficulty in applying the method would be encountered with regard to 
satellite measurements because of their relatively low resolution. Instruments 
such as the AVHRR ( 1 km resolution ) may be incapable of producing a wide 
range of NDVI and surface temperature over a target area. Thus, it may be 
necessary to extrapolate the axis of variation and to composite satellite mea- 
surements over successive images. Nemani and Running (1989b) have re- 
cently presented some encouraging results showing that the AVHRR is capa- 
ble of determining an axis of variation over forests. Their figures show that 
the slope of the axis of variation varies according to the amount of surface 
moisture. 

SUMMARY 

This study suggests a method to derive a spatial distribution of fractional 
vegetation cover and soil moisture in the surface and root zone over patchy 
or sparse vegetation. The method makes use of two types of image products: 
the distribution of NDVI vs. radiometric surface temperature (the NDVI/ 
T c ) and the “arch” diagram, which shows the distribution of radiometric sur- 
face temperature in subgroups of pixels vs. the standard deviation of the ra- 
diometric surface temperature in those subgroups. These two diagrams aid in 
identifying asymptotic limits of the sunlit leaf and the sunlit bare soil temper- 
ature. It also allows one qualitatively to assess the degree of noise produced 
by small-scale variations in soil moisture and leaf shading. 

The boundary layer model is used to estimate soil surface and root zone 
water contents, given the asymptotic vegetation and bare soil temperatures. 
A crucial aspect of the method is to define the relationship between NDVI 
and surface radiometric temperature, which we call axis of variation. From 
this, we obtain average soil moisture values for surface and root zone, and a 
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vegetation fraction, surface turbulent heat fluxes and LAI as functions of 
NDVI for each pixel. 

Some subjectivity is inherent in the method. An axis of variation may not 
always be easy to define, for example when small-scale variations in soil mois- 
ture are large (perhaps owing to irrigation), where the amount of bare soil is 
negligible or where the soil surface temperature and leaf temperatures are 
similar, as in the case of wet soil. Large horizontal variations in soil moisture 
increase the scatter. Nevertheless, even when the axis of distribution is ob- 
scured because of noise, the NDVI/r c and arch diagrams yield useful quali- 
tative information on the distribution of vegetation cover, small-scale varia- 
tions in soil moisture and vertical gradient of soil water content. A reliable 
axis of variation could be determined by compositing measurements ob- 
tained from a sequence of images. 

The method needs to be tested on many types of soils and vegetation. Ap- 
plication to analysis of AVHRR satellite data may prove useful over different 
types and mixes of vegetation, although it is still too early to determine if 
lower-resolution satellite radiometers can define a usable axis of variation. 
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ABSTRACT 


TRANSIENT WATER STRESS IN A VEGETATION CANOPY: 
SIMULATIONS AND MEASUREMENTS 

Measurements of evapotranspiration and radiometric surface temperature, 
supplemented by model simulations, illustrate the effect of transient water stress on the 
evapotranspiration and radiometric surface temperatures over a corn ( Zea Mays L. ) 
canopy. Transient stress is manifested by a plateau in the evapotranspiration and a 
warmer radiometric temperature of the canopy. The cause of the evapotranspiration 
plateau is not solely due to reduced soil water content but to excessive atmospheric 
demand in the face of internal plant constraints. One of these constraints is a threshold 
epidermal leaf water potential, below which stomatal resistance increases very rapidly with 
decreasing leaf water potential. Accordingly the evapotranspiration plateau can occur at 
relatively high soil water content. 

The sensitivity of the canopy temperature to changing soil water content is weak at 
high soil water content. The evapotranspiration plateau and its effect on radiometric 
surface temperatures would not be detectable until, in this case, the soil water content 
decreases to about half of field capacity. Since evapotranspiration collapse (wilting) 
would occur at about 35% of field capacity, it may be possible to remotely monitor the 
onset of crop water stress over dense vegetation before catastrophic damage occurs. 
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Introduction 


A problem in remote sensing is that the infrared temperature of dense vegetation is 
rather insensitive to soil water content. A current ("root-shoot") controversy amongst 
plant scientists (Kramer, 1988; Passioura, 1988; Boyer, 1989) may have important 
implications for solution of this problem. 

Let us very briefly summarize the two positions. Adherents of root control cite 
laboratory evidence showing that stomatal resistance varies according to signals carried in 
sap from roots to leaves by a hormone (ABA), called abscisic acid. These root signals 
inform the leaves of the soil water status and allow the stomates to adjust accordingly 
(Gollan et al, 1986; Schulze et al., 1987; Munns and King, 1988). 

A counter - argument in favor of leaf control of stomatal resistance offers evidence 
that the epidermal cells in the leaf are able to sense water stress. Leaf water stress is more 
closely related to the leaf turgor and hydraulics of water flow to the leaf than to 
biochemical signals from the roots (Frensch and Schulze, 1988; Kramer, 1988). 
Resolution of this issue awaits more research. 

Alternately stated, the control of stomatal resistance by the leaf suggests that the 
former is related to leaf water potential. This implies that the leaf can react to transient 
periods of water stress brought about by high atmospheric demand and a loss of positive 
turgor pressure in the leaves (Jones, 1983). Transient periods of water stress could then 
occur in principle even when the soil is wet, provided that the atmospheric demand 
exceeds the plant's capability to supply water from its roots. In such a case, transient leaf 
water stress would be manifested by a midday plateau or dip in the transpiration, with an 
associated rise in leaf temperature and sensible heat. 
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The onset of such transient leaf water stress can be modeled by relating epidermal 
leaf water potential to stomatal resistance and invoking an epidermal leaf water potential 
threshold ( \p c ), below which stomatal resistance increases rapidly with decreasing leaf 
water potential. 

The concept a leaf water potential threshold originated with Cowan (1965), who 
defined it in terms of a supply function. Turner (1970; 1974) and many others (e.g. Nizinski 
and Saugier, 1989) have presented striking observational evidence that such a leaf water 
potential threshold occurs in a variety of crops. The leaf water potential threshold 
constitutes an important parameter in the plant component of some boundary layer or 
climate models (Wetzel and Chang, 1987; Sellers and Dorman, 1987; Lynn and Carlson, 
1990). 


The purpose of this paper is to [1] present additional observational and modeling 
evidence of transient water stress, [2] illustrate the effects of the transpiration plateau on 
the canopy radiometric temperature and [3] describe the factors responsible for the onset 
of the transpiration plateau, such as soil moisture. This paper also considers the point at 
which the transient stress can be detected by remote measurement of surface 
temperature, including the effects of surrounding bare soil, cloud cover and how these can 
mask the stress signal. 

2. Field Measurements 

Micrometeorological field measurements were made over a two - year period 
(1987 - 88) at Rock Springs, a research agricultural site near The Pennsylvania State 
University. These measurements consist of net radiation, air temperature, wind speed and 
relative humidity, taken from a tower at approximately 9 meters above the ground. 
Surface radiometric temperature was measured by a radiometer (Optitherm II) mounted 
on a second tower and trained towards the surface at an elevation of about 3 meters at an 
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angle of 20 degrees from the vertical. Surface sensible heat flux was obtained using two 
methods: [1] Directly from the sonic anemometer and [2] from bulk aerodynamic 
formulae using measurements of surface radiometric and air temperatures, wind speed 
and crop height. Surface roughness and displacement height was obtained using standard 
formulae. Ground heat flux was determined in two ways: A flux plate embedded at a 10 
centimeter depth and using temperature profiles obtained from ground thermistors buried 
up to a depth of 40 centimeters. Soil water content was measured using gypsum blocks 
buried at 5 substrate levels from 2.5 to 40 centimeters in depth. Instruments operated 
continually and the measurements were recorded as half - hour averages on computer 
diskettes. 

Surface evapotranspiration was determined from the surface energy budget by 
subtracting the sum of the ground flux and the surface sensible heat flux from the net 
radiation flux. During the two growing seasons, there was only one period in which, [1] the 
measurements functioned continuously for several days under sunny skies, [2] the 
vegetation was dense enough to permit the assumption that evaporation flux from the 
underlying soil was relatively small compared to the total evapotranspiration and [3] the 
sensible heat fluxes measured independently by two methods were in very close 
agreement. This period was during July 1987 in a corn canopy ( Zea Mays L. ), which was 
about 2.5 meters high and had a leaf area index of about 7 (Belles, 1990). 

These measurements nonetheless provide good examples of the evapotranspiration 
plateau. Figure 1 shows this plateau, which is indicated by a sudden change in slope of the 
evapotranspiration curve. The plateau lengthened progressively between July 18^\ when 
it lasted about two hours, and July 22 n< ^, when it lasted from about 0900 until 1500 local 
time. 


Figure 2 shows the progressive decrease in the evapotranspiration during the 



middle part of the day from July 18 th to 22 th . Differences are slight during the morning 
and late afternoon before and after the time of the plateau. Measured differences from 
day to day were also largest near midday. Similarly, the radiometric minus air 
temperature differences (Tj r - T a ) between July 18^ and July 22 n< ^ were largest (about 
3°C) near midday and small during the early morning and late afternoon (Fig. 3). 

The increasingly early onset and lengthening duration of the evapotranspiration 
plateau during these four days is due to a progressive drying of the soil. Indeed, soil water 
contents (Fig. 4) at both 2.5 and 20 centimeters diminished from values in excess of 0.2 (by 
volume) to less than 0.1 between July 4^ and July 30^, with the surface drying out at a 
greater rate than that at the deeper level. More precisely, the deeper layer, which we 
associate with the root zone layer, dried out from 0.74 of field capacity (0.34 by volume 
for the soil type - Hagerstown loam) on July 4th to 0.5 of field capacity on July 18th to 0.28 
of field capacity on July 30th. 

3. Simulations - boundary layer model with plant component 

A tool used to interpret our measurements is a boundary layer model with plant 
component. The model is able to simulate the evapotranspiration plateau from which we 
can analyze factors responsible for onset and duration of the plateau. The atmospheric 
boundary layer model has been described by Carlson (1986), its vegetation component by 
Taconet et al. (1986) and the plant hydraulics by Lynn and Carlson (1990) and Carlson 
and Lynn (1990). We will restate in this paper, only those aspects of the vegetation model 
that pertain to the simulation of the plateau. 

The model calculates evapotranspiration from a balance between the flow of water 
from the substrate to the plant (F w ), the transpiration (vapor) flux from the leaves to the 
surrounding atmosphere (L e Ef) and the water flux from storage (q s ). 
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The equation is expressed as: 

L e E f = F w + % Ha] 

For steady state (q s = 0) this water balance can be expressed in resistance notation as: 


~ 'h, ~ H P c p 
z g + z p " 7 

Here, F w is the flux of liquid water from the soil to the leaf. It is expressed by a drop in 
water potential from the soil around the roots (^g) to the leaves in the mesophyll (i£j) 
divided by the sum of the soil plus plant resistances ( Z g + Z p). This is equated with 
transpiration (L g Ef), the flux of water vapor (latent heat) from the leaves to the air. It is 
given by the vapor pressure drop from the leaf ( ej at the temperature of the leaf Tj ) to 
the surrounding air (e a f) divided by the sum of leaf plus its boundary layer resistance 
( q + r a f ). H is the gravitational potential ( H = p w gh, where p w is the liquid water 
density, g is the gravitational constant and h represents the height of the leaves above the 
roots ). The psychrometric constant is 7 ( 7 = PCp/(0.622L e ), where P is the atmospheric 
pressure, Cp is the specific heat at constant pressure and L g is the latent heat of 
condensation. 

Leaf resistance consists of stomatal and cuticle resistances, which exist in parallel. 
Since we set the cuticle resistance to be much larger than the stomatal resistance, we will 
refer to the leaf resistance rj as the average stomatal resistance r s ; the latter is defined for 
a sunlit leaf near the top of the canopy. 

Substrate water is accounted for in two layers, a surface layer and a root zone layer. 
The former controls the evaporation from the soil surface beneath the plant and the latter 
governs i£g and Z g. In this case since there was dense vegetation cover we made the 
assumption that evapotranspiration and transpiration are virtually identical quantities. 
The ground water potential (i^g) and the soil resistance ( Z g) are both dependent on the 


e i(Tj) - e^ 


r s + r af 


[lb] 
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soil water content and soil type. We calculate the soil water potential and soil resistance 
using simple formulae given by Cosby et al. (1984). The evaporative flux is computed for 
the underlying soil layer. 

The stomatal resistance (r s ) is assumed to depend on leaf water potential, vapor 
pressure deficit and sunlight intensity. We use the following relationship to relate these 
parameters to stomatal resistance: 

r s - r min * £(* e ) f < s > P) 

where f(S) and f(^ e ) are, respectively, functions of incident solar flux and epidermal leaf 
water potential and r^Q is a minimum stomatal resistance. The two functions, f( \p e ) and 
f(S), vary from one to infinity; is the value of r s when the functions are both equal to 
one. At f(S) = 1.0, the solar flux is assumed to be a saturation value (S 0 ), which we 
assume to be 1000 Wm'^. At f(^ e ) = 1.0, the leaf is saturated with water and \p Q = 0. 

The effect of vapor pressure deficit on stomatal resistance is accounted for by 
setting a proportionality between epidermal minus mesophyll water potential difference 
and the vapor pressure deficit (e^Tj) - e af ). Lynn and Carlson (1990) show that this 
formulation realistically captures feedbacks between vapor pressure deficit and stomatal 
resistance. However, since the stomates in corn are thought to possess a low sensitivity to 
vapor pressure deficit (B. Choudhury, private communication), we let ^ 

The functions f(iA e ) and f(S) are highly exponential, first increasing slowly from a 
value of one and then much more rapidly with decreasing values of \(> e or S. To capture 
this non-linearity in the simplest possible manner consistent with measurements, we 
separate each exponential function into to two linear regimes, one with a small slope and 
the other with a large slope. Such a model, called 'linear discontinuous' by Fisher et al. 
(1981), has also been adopted by Nizinski and Saugier (1989). For the leaf water potential 
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function, we write: 


f(*e> ='W r min < r min + Ve + b 2^e * *c) ** > ! 3 1 

where and b£ are constants representing the respective slopes of the two lines and \p c is 
the ’threshold leaf water potential’, which lies at the intersection of the two straight lines. 
The operator 5 ^ is equal to zero when is in the sub - threshold (slowly changing) regime, 
(larger than \p c ). When reaches \p c or falls below it, 6^, = 1.0 and <^ e equals <P C in the 
term multiplied by bj. 

We found that the best fit with measurements was to set bj equal to zero for all 
simulations (Table 1). In so doing, we imply that there is no sensitivity of stomatal 
resistance to leaf water potential until the threshold is reached. Thus, when the leaf water 
potential is above threshold, substrate water content affects the water balance equation 
(Eqn. lb) only by way of the ground water potential ig and the ground resistance Zg, both 
of which change very slowly with time. When the leaf water potential decreases below \p c , 
the contribution by the \>2 term is very large (Table 1) and the stomatal resistance 
increases rapidly with decreasing leaf water potential. 

The threshold leaf water potential accompanied by a large magnitude of b 2 
produces a transpiration plateau. The dependence of r s on leaf water potential appears to 
be necessary for simulating this effect. This is the case even when stomatal resistance is 
linked to vapor pressure deficit. 

A similar linear discontinuous function is used for f(S). Since the threshold solar 
flux (S c ) is usually rather low ( 50 Wm* 2 ; Table 1 ), the behavior of f(S) is unimportant for 
this study and will not be discussed further. 
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Water storage also can affect the shape of the transpiration plateau (Carlson and 
Lynn, 1990). The role of capacitance is relatively minor however, and thus we will not 
discuss the nature of the flux q s term (Eqn. la) except in passing. Capacitance is 
calculated in the model but it is omitted in (Eqn. la) because a discussion of its 
mathematics is felt to be unnecessary for this paper. 

Equations [1] - [3] constitute a quadratic which can be solved analytically for the 
leaf water potential ( Lynn and Carlson; 1990 ) and by substitution, the transpiration. 
These equations form part of a larger boundary layer model. The larger model operates 
simultaneously in two modes, one for vegetation with underlying bare soil and the other 
for bare soil without vegetation. Bare soil and vegetation components are blended to 
account for situations of partial vegetation cover ( Carlson et al., 1990 ). For the most 
part, we assume full vegetation cover and suppress the separate bare soil component; 
however some simulations were made for partial vegetation cover with a specified 
vegetation fraction. 

Calculations of r a f and e a f, Tj and solar flux were made using equations discussed 
by Taconet et al. (1986). Fluxes computed in the model pertain to a single leaf. Extension 
of these fluxes from a leaf to a canopy is currently undergoing considerable scrutiny 
(Baldocchi et al., 1990; Rochette et al., 1990; Mascart et al., 1990). Our method for 
scaling from a leaf to a canopy is to multiply the individual leaf fluxes by a shelter factor to 
leaf area index ratio. The shelter factor accounts for the fact that not all leaves transpire 
at the sunlit amount due to fact that solar radiation decreases with height beneath the top 
of the canopy. For dense vegetation, the shelter factor to leaf area index ratio approaches 
two. The reader is referred to Taconet et al. (1986) and Lynn and Carlson (1990) for a 
more detailed discussion of our scaling factor. 
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4. The transpiration plateau 

Simulations of evapotranspiration were made for July 18 1 * 1 , 21 st and 22 nc *, 1987 
(Fig. 1) using values of the vegetation parameters shown in Table 1; these are similar to 
values adopted by Lynn and Carlson (1990) for their corn simulations. Soil water contents 
for the surface and root zone are given in Figure 4. Except for the soil water contents and 
the initial atmospheric conditions, the model parameters are identical for all three days. 

The purpose in performing these simulations is to demonstrate that the model is 
capable of reproducing the transpiration plateau during the period 18^ to 22 n ^ July. We 
obtained close agreement with the measurements for values of threshold leaf water 
potential and other parameters in Table 1. The threshold water potential of \p c = -15 bars 
is typical for temperate latitude crops (Nizinski and Saugier, 1989). Agreement between 
measurement and model results depends upon the choice of Zp, which was 0.06 bars 
(Wm'2)*l ( Table 1 ). A smaller value of Zp would have required our using a smaller 
value of the threshold leaf water potential. Lynn and Carlson (1990) determined a 
threshold leaf water potential of -13 bars and Zp = 0.04 bars (Wm'^)^ for com. 

The evapotranspiration plateau tends to slope downward with time. We can 
simulate this slope in two ways: [1] By the inclusion of plant capacitance and [2] by 
relating stomatal resistance to vapor pressure deficit. Since the vapor pressure deficit 
effect was suppressed in the simulations, the slope in the modeled plateau was due to 
capacitance. 

A primary objective of this paper is to show how the plateau would have evolved 
during a period of soil drying, such as between the 4 1 * 1 and 30^ of July, 1987. To 
accomplish this, we used the soil water contents obtained from inspection of Figure 4 for 
seven days during the drying period ( July 4^, 10^, 18^, 21 st , 22 n< *, 26^ and 30^ ). We 
further assume that these days were clear, which was generally not the case except for July 
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18 th , 21 st and 22 nd . Initial atmospheric conditions for these three days were obtained 
from the Pittsburgh upper air sounding ( 200 km from the measurement site ) and from 
local surface weather observations. 

Results proved to be insensitive to initial atmospheric conditions but highly 
sensitive to initial deep layer soil water content. In view of the uncertainty in specifying 
initial atmospheric conditions (due to the absence of local rawinsonde measurements), we 
used the initial atmospheric conditions for July 22 nd in all simulations except for July 18^ 
and 21 st , when field measurements existed for testing the model. Our intention is to 
investigate the transpiration plateau, showing its evolution on clear days in response to a 
progressive drying of the soil. 

Figure 5 is an example of a simulated decrease in the evapotranspiration as a 
function of soil water content, which we specify as a function of surface (f) and deep soil 
(fsub) water contents, expressed as fractions of the field capacity. Note that the two 
wettest cases are almost identical. Lynn and Carlson (1990) also demonstrated that very 
little change occurs in the transpiration over corn until the soil water content decreases to 
half of field capacity. In these simulations an obvious plateau occurs when the root zone 
soil water content is decreased below 59% of field capacity, which appears after the 18 th . 
A collapse in the transpiration happens when the ground water potential ^ is reduced to 
the threshold leaf water potential (-15 bars), which is equivalent to a soil water content 
of 0.12 by volume ( fsub = 0.35 ). 

Evapotranspiration for the two wettest days in Figure 5 also shows a slight 
flattening near midday as the leaf water potential briefly reaches the threshold potential. 
Figure 6 suggests that a small plateau was present some time prior to the 18 , possibly as 
early as July 4 1 * 1 , when the soil water content was 74% of field capacity. Thus, the 
threshold potential may be reached in relatively well - watered vegetation providing the 
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atmospheric demand is sufficiently large. 


An important aspect of both Figures 2 and 5, which has also been pointed out by 
Lynn and Carlson (1990), is that the variation in transpiration and radiometric surface 
temperature of the corn crop canopy is small from day to day when there is no plateau. 
Once the plateau forms, changes in evapotranspiration and radiometric temperature from 
one day to the next increase during the time when the plateau is present. By setting bj to 
a value of zero in Eqn. 3 the evapotranspiration curves are virtually identical until the 
onset of the plateau. Stated alternately, the curves of Figure 5 are all alike before the 
onset of the plateau because leaf water potential is decoupled from stomatal resistance in 
the sub - threshold regime ( ). 

5. Radiometric temperature of the canopy 

Radiometric surface temperatures reflect the presence of the transpiration plateau. 
To remove the variable effects of the atmosphere on radiometric surface temperature, we 
compute the radiometric surface and ambient air temperature difference. Simulations 
(Fig. 7) and measurements (Fig. 3) of the T if - T air differences increase most rapidly with 
decreasing soil water content during the period of the transpiration plateau. We are 
unable to explain the reason for the higher and more rapid increase in the measured 
versus simulated surface minus air temperature differences during this four - day period. 
While the measured values seem somewhat large, being 1 - 2 °C positive at radiation 
sunset and sunrise (Fig. 3), it should be noted that measured sensible heat fluxes for both 
the aerodynamic method and sonic anemometer were in close agreement. The important 
point, however, are that measurements and simulations exhibit similar response to soil 
drying, particularly with regard to the midday increase of radiometric surface temperature. 

This response of leaf temperature to threshold leaf water potential is evident in 
simulated Tj r - T a j r throughout the entire drying period (Fig. 8). Largest increases from 
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day to day occur after the plateau has formed. The two wettest cases, corresponding to 
4 1 ^ and 10 l “ July, show virtually no differences at any time during the day. In the driest 
case, Tj r - T a j r is at a maximum and corresponds to the evapotranspiration collapse 
( re: Fig. 5 ). 

As the length of the plateau increases, so does the midday radiometric surface 
temperature. Figure 9a demonstrates a radiometric temperature increase almost linearly 
with plateau length. The increase in radiometric surface temperature following July 4 1 * 1 
remains less than a degree during the first two weeks. When the plateau length exceeds 
about 200 minutes, the increase in midday radiometric surface temperature with 
decreasing soil water content becomes noticeably more rapid with time. Similar behavior 
is exhibited in the measurements (Fig. 9b), although Tj r - T a j r is larger and the plateau 
length somewhat shorter than for the simulations. 

Simulated and measured T ir - T a ^ versus soil water content are given in Figure 10. 
Again, the measured radiometric surface temperatures are higher than those simulated, 
but more importantly, the midday surface radiometric temperature first increases 
gradually and then considerably more rapidly after July 18^. If one concedes that a 
reasonable measurement uncertainty for radiometric measurement of a vegetation canopy 
is at least plus or minus 1°C, Figures 9 and 10 suggest that one is unable to detect a 
reduction in soil water content below field capacity for a period of at least two weeks after 
July 4*k, by which time the soil water content diminished to about half of its field capacity. 

6. Effects of vegetation cover soil and cloud cover 

(a) Vegetation cover 

The effect of bare soil patches is to mask the evapotranspiration plateau. We made 
simulations in which fractional vegetation cover was varied, while maintaining the same 
surface and root zone water content values obtained from Figure 4 ( Figs. 11a & b ). 


12 



Increasing the fraction of exposed bare soil causes the plateau to decrease and 
subsequently vanish when the vegetation cover is reduced to 50% of full cover (Fig. 11a). 
At 50% vegetation cover, Tj r - T a j r is larger and exhibits a steeper increase with 
decreasing soil water content than for full vegetation cover (Fig. lib). Although these 
changes in temperature for partial vegetation cover are more easily detectable by a 
radiometer, they represent a mixed signal from both surface and root zone. 

(b) Cloud cover 

In addition, cloud cover also affects the transpiration plateau. A reduction of solar 
flux due to cloud cover does not produce a corresponding reduction in the transpiration 
flux. Decreasing the solar flux by as little as 25 %, simulating scattered cloud cover, 
effectively eliminates the plateau (Fig. 12) and the corresponding rise in T ir - T air 
(Fig. 13). The disappearance of the water stress signal occurs as a result of a reduction in 
the solar flux and therefore in the atmospheric demand. Katerji and Hallaire (1984) also 
show that a reduction in solar demand reduces the decrease in the leaf water potential. A 
further decrease in soil water content, however, would cause the plateau to reappear for 
scattered and eventually for broken sky cover, provided that is sufficiently low. 

7. Conclusions 

Although soil moisture limitation constitutes an underlying cause of the 
evapotranspiration plateau, its onset, duration and shape are determined by a combination 
of factors, particularly those of atmospheric demand and internal plant constraints. 
Demand constitutes the incident solar flux and the vapor pressure deficit. The important 
plant constraint is the threshold leaf water potential, which is necessary for simulating the 
evapotranspiration plateau. Given that stomatal resistance increases rapidly for 
decreasing leaf water potential below the threshold, simulations with the plant model 
indicate an increase in the intensity of solar radiation, a decrease in the ambient relative 
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humidity or a decrease in soil moisture contribute to an earlier formation and longer 
duration of the evapotranspiration plateau. Therefore, the evapotranspiration plateau 
may occur when the soil is relatively moist, providing the atmospheric and solar demand 
are high. It may not occur when the soil is dry if the demand is weak. 

Onset of the evapotranspiration plateau is manifested by an increase in the leaf 
( radiometric surface ) temperature and sensible heat flux. Our simulations suggest that 
the largest changes in radiometric surface temperatures occur during the period of the 
plateau, notably near 1300 local time. We show that an increase in radiometric canopy 
temperature with decreasing soil water content below field capacity is detectable by a 
radiometer only when the plateau lasts for more than a couple of hours. In the case of a 
dense corn canopy, detectability begins when the soil water content decreases to about 
50% of its field capacity. Since the true wilting point for the corn was about 35% of field 
capacity ( essentially the point at which the soil water potential reaches the threshold leaf 
water potential ), remotely sensing the onset of severe crop water stress may be possible 
using radiometric surface temperatures. 

Finally, we surmise that the evapotranspiration plateau operates to protect the 
plant against excessive water loss. However, the plateau implies an increase in the leaf 
temperature, likewise in the vapor pressure deficit between leaf and air, and in stomatal 
resistance. This would lead to a further destructive escalation in leaf temperature. In view 
of the adverse effects of such a positive feedback, it seems reasonable to suppose that the 
evapotranspiration plateau may not constitute a routine and optimum protective strategy 
for a plant under normal conditions. Rather, we feel that plants have evolved such a 
survival mechanism in response to their indigenous habitat such that they experienced the 
evapotranspiration plateau only under conditions of abnormal soil or atmospheric 
dryness. 
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Table 1 


List of some parameters used in the simulations 


Parameter 

Value 

Units 

^min 

25 

s nr 1 

^c 

-15 

bars 

bi 

0 

s m* 1 bars' 1 

b 2 

500 

s m' 1 bars' 1 

S c 

50 

s m* 1 W ' 1 m 2 

Z p 

0.06 

bars (Wm* 2 ) - 1 
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The measured latent heat (evapotranspiration) fluxes (Wm'^j squares) 
from aerodynamic method (over com; leaf area index = 7.19) based on 
the surface radiometric temperatures at Rock Springs, 18, 21 st and 
22 nd July, 1987 and the simulated fluxes (triangles) for the same days. 
Vertical arrows denote duration of the evapotranspiration plateau in the 
measurements. 

Measured evapotranspiration fluxes (Wm'^) versus time for 18 , 21 st 
and 22 nd July, 1987. 

Measured surface radiometric minus air temperature differences (°C) 
versus time for 18^, 21 st and 22 nd July, 1987. 

Volumetric soil water content at 0700 local time during July 1987 as 
measured by gypsum blocks at depths of 2.5 cm (squares) and at 20 cm 
(triangles). The dashed horizontal line denotes the volumetric field 
capacity for the soil and the hanging bars are proportional to the daily 
rainfall, which is measured against the scale in centimeters at the right. 
Vertical arrows along horizontal scale indicate dates in July as referred 
to in the text. 

Simulated evapotranspiration fluxes (Wm ) versus time but for 
differing surface and root zone soil water concentrations obtained from 
Figure 4 (legend shows surface (f) and root zone (fsub) water contents 
expressed as fractions of the field capacity, which is 0.34 by volume). 
From top to bottom these curves represent hypothetical clear sky days: 
July 4*k, 10 th , 22 nd , 26^ and 30^. ( All initial atmospheric conditions 
pertain to July 22 nd . ) 
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Volumetric root zone water content ( fsub; fraction of field capacity ) 
versus length of the evapotranspiration plateau (min) for simulations of 
hypothetical clear sky days ( squares ) and for measurements ( crosses ). 
Numbers aside data points denote the dates in July 1987. 

Simulated surface radiometric minus air temperature differences (°C) 
versus time for three days in July 1987 ( see legend ) for soil water 
contents obtained from Figure 4. 

Simulated radiometric surface minus air temperatures differences (°C) 
with time for the five idealized clear sky days shown in Figure 5. 

Surface radiometric air temperature differences (°C) at 1300 local time 
versus length of the evapotranspiration plateau (min) for simulations 
(solid line and squares). Numbers refer to dates of hypothetical clear 
sky days in July 1987. 

Same as Figure 9a but a comparison of measured (dashed line and 
crosses) and simulated (solid line and squares) evapotranspiration 
plateau lengths. 

Radiometric surface minus air temperature differences (°C) at 1300 
local time versus root zone volumetric water content (cnr/cm J ) for 
simulations (lower curve) and measurements (upper curve). 

Simulated variation of evapotranspiration fluxes (Wm'^) versus time for 
differing fractional vegetation cover. 
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Figure 13. 


Radiometric surface minus air temperature differences (°C) at 1300 
local time versus root zone volumetric water content (cm J / cm J ) for full 
vegetation cover (lower curve) and 50% vegetation cover (top curve). 
Numbers refer to dates of idealized clear sky days in July 1987. 

Radiometric surface minus air temperature differences (°C) versus 
fractional root zone water content (fsub) simulated for 100% and 50% 
fractional vegetation cover using initial atmospheric conditions of 22 n( ^ 
July 1987. 

Simulated evapotranspiration fluxes (Wm'^) versus time for clear, 
scattered and broken sky cover; these conditions correspond, 
respectively, to full sun, 25% and 40% reduction in solar flux using 
initial atmospheric conditions of 22 n< * July 1987. 
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Figure 2. Measured evapotranspiration fluxes (Wm'^) versus time for 18^, 21 st and 
22 nc * July, 1987. 
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Figure 3. Measured surface radiometric minus air temperature differences (°C) 
versus time for 18 th , 21 st and 22 nd July, 1987. 


Tir-Ta vs Plateau Length 

July 1987 



Plateau Length (min) 


Figure 9b. 


Same as Figure 9a but a comparison of measured (dashed line and 
crosses) and simulated (solid line and squares) evapotranspiration 
plateau lengths. 
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Figure lib. 


.adiometric surface minus air temperature difference^(°Q at 1300 
>cal time versus root zone volumetnc water content (cm /cm ) tor tun 
egetation cover (lower curve) and 50 % vegetation cover (top curve), 
lumbers refer to dates of idealized clear sky days m July 1987. 
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Figure 4. Volumetric soil water content at 0700 local time during July 1987 as 
measured by gypsum blocks at depths of 2.5 cm (squares) and at 20 cm 
(triangles). The dashed horizontal line denotes the volumetric field 
capacity for the soil and the hanging bars are proportional to the daily 
rainfall, which is measured against the scale in centimeters at the right. 
Vertical arrows along horizontal scale indicate dates in July as referred to 
in the text. 
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Figure 5. 


Simulated evapotranspiration fluxes (Wm*^) versus time but for differing 
lurface and root zone soil water concentrations obtained from Figure 4 
legend shows surface (f) and root zone (fsub) water contents expressed as 
Tactions of the field capacity, which is 0.34 by volume). From, top to 
>ottom these curves represent hypothetical clear sky days: July 4 , 10 l , 
nna 26 th and 30 th . ( All initial atmospheric conditions pertain to July 


>nd’ 
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Figure 6. Volumetric root zone water content ( fsub; fraction of field capacity ) 

^ versus length of the evapotranspiration plateau (mm) for simulations of 

hypothetical clear sky days ( squares ) and for measurements ( crosses ). 
Numbers aside data points denote the dates in July lyo/. 
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Figure 7. 


imnlated surface radiometric minus air temperature differences (°C) 
S toe fonteee C in July 1987 ( see legend ) for soil water contents 


obtained from Figure 4. 
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Figure 8. Simulated radiometric surface minus air temperatures differences ( °C) with 
c ” time for the five idealized clear sky days shown in Figure 5. 
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Surface radiometric air temperature differences (°C) at 1300 local time 
versus length of the evapotranspiration plateau (min) for simulations 
(solid line and squares). Numbers refer to dates of hypothetical clear sky 
days in July 1987. 
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Figure 10. Radiometric surface minus air temperature differences (°C) at 1300 local 
& time versus root zone volumetric water content (cm /cm ) tor 

simulations (lower curve) and measurements (upper curve). 
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Figure 11a. Simulated variation of evapotranspiration fluxes (Wm ) versus time for 
differing fractional vegetation cover. 
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Figure 13. Simulated evapotranspiration fluxes (Wm"^) versus time for clear, 
scattered and broken sky cover; these conditions correspond, respectively, 
to full sun, 25% and 40% reduction in solar flux using initial atmosphenc 
conditions of 22 nd July 1987. 






